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INTRODUCTION 



I . 



A. BACKGROUND 

The stream of data generated during an anti-air warfare 
(AAW) operation at sea can be overwhelming, and an efficient 
and complete analysis of the information provided by all the 
different shipboard sensors will take more time than 
available. One solution to the problem is the design of a 
computer based Decision Support System ( DS3 ) . This system will 
have the capabilities of today's computer technology, 
especially graphics, so that it will be able to utilize 
information gathered on own ship, on enemy forces, and also 
historical and intelligence data. The system will effectively 
predict enemy near term moves, will organize all of this 
collected information, and will present it in a form 
(pictures) that will be easy to use and understand. With a 
sophisticated algorithm, it can lead to appropriate decision 
making. 

One objective of this thesis is to design a small 
interactive DSS and to show examples of a few modules written 
to deal with various features of the engagement of an 
incoming threat. After analyzing the available information 
about this target and about our weapons, an optimized use of 
our power against this target will be suggested to the 
user. An analysis of his action will be carried out and 
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the results displayed. Another objective is to show examples 
of how to convert some of the raw information provided by 
some of the ship borne sensors to a comprehensive picture 
on the screen of an IBM PC/AT computer. The scope of this 
thesis will be limited to an AAW operation conducted from the 
Combat Information Center (CIC). This work does not provide 
an operational capability because it is meant to be 
educational. Most of the elements used here are generic since 
the real ones are classified. 

B. I3EAS 

The Interactive Simulation of an Engagement At Sea (ISEAS) 
is a PC based computer program, under devel opement at the 
Naval Postgraduate School , where many of the features of a 
real engagement are simulated. It is written in the C language 
and reproduces the functions of a CIC, from the use of the 
radars to the direction of the weapons. 

This paper has two main parts. The first part (Chapters. 
II, III and IV) deals with the tactical situation and the 
decision-making. Thus, a CIC doctrine is simulated given the 
information provided by the RADAR module and the WEAPONS 
module which are written by other students. The second part 
(Chapters V and VI) will look at some of the ways to display 
information with the aid of graphics. 
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II. THE BIG PICTURE 



A. DEFINITIONS 

What happens when there is a target to engage? How is the 
engagement decided, and when? What happens during the 
engagement, and afterwards? What is the role of the Combat 
Information Center? How does it conduct this task? This 
chapter will introduce the user of ISEAS to the different 
events that occur during the engagement and to the different, 
possible outcomes. Also, the user will be made familiar with 
the structure of ISEAS and be able to identify the different 
modules and understand their relationship to each other. 

A big picture of the engagement will be drawn here giving 
a clear idea of how an outcome can be reached. To do so, an 
Essential Elements Analysis is conducted, and a relationship 
between a wanted outcome and an end-event is established. 

First, the definitions of some keywords used throughout 
this thesis are needed. Usually, these keywords mean different 
things because of the variety of the contexts and also because 
the tactical doctrines keep changing with time. The main 
difference is their scope and extent of use. The author was 
unable to find any specific written definitions. But here they 
are defined and the mean the following. 

Combat Information Center (CIC): it is the space within 
the ship in which the battle is fought. Tactical control of 
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all the ship's combat systems resides in the CIC. Inter-ship 
communications and other command and control functions are 
also performed in the CIC. Specifically, the following five 
functions are performed: 

1. Information collection 

2. Information processing 

3. Information display 

4. Information dissemination 

5. Weapons direction (involving the fire control 

system) 

Fire Control System ( FCS ) : it is the set of computers, 

consoles, data buses interfacing with sensors, weapons 
launchers and weapons direction equipment. It takes inputs 
from the sensors, processes them in the computers and 
generates firing solution computations and weapons direction. 

Command and Control unit: it is the processor of 

information that generates solutions other than the firing 
computations, such as information about a detected target, 
pulled out from a large database system; or a suggestion of an 
optimised use of the ship's weapons against a raid. Such 
suggestions are based on a set of rules of different tactical 
situations and doctrines, implemented in the algorithms used 
by the processor. 
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B. ESSENTIAL ELEMENTS ANALYSIS 



In a one-on-one engagement between a ship and an aircraft 
or a missile there are four possible outcomes: 

1. The target is killed, and the ship survives 

2, The target survives, and the ship is killed 

3. The target survives, and the ship also survives 

4, Both the target and the ship are killed 

An Essential Elements Analysis ( EEA ) can be conducted for 
each outcome to see what are the events that lead to it. 
However, since we are more interested in looking at the 
action of the ship against the attacker, and since we assume 
(at the present time) that our ship is invulnerable, an EEA 
for the target kill alone will be conducted. 

There are two possible outcomes for the target: 

1 . Target ki 1 led 

2. Target not killed 

A target is killed if the ship succeeds in countering it 
and denying its penetration at one of the defense envelopes 
(long range, medium range or short range missile envelope); or 
if the target succeeds in entering all the envelopes but fails 
to function accurately and misses the ship. Again since our 
concern is to actual ly' prevent a successful enemy attack, we 
will look at the first possibi lity; the ship effectively 
counters the target. For the ship to succeed in killing the 
target the following events should occur successfully: 

- target detected in time 
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- target identified and designated correctly as enemy 

- target assigned to the fire control system 

- target acquired by the fire control system 

- target tracked accurately 

- a fire control solution generated correctly 

- a proper weapon is available and ready for use 

- missile is launched 

- missile intercepts the target 

- target is killed by missi le' warhead 

The algorithm of the engagement of a target by a weapon 
will follow these straightforward steps, provided that each 
step is accomplished. An assessed probability of happening can 
be computed or determined from a table prior to the actual 
engagement. The ultimate result is a value of how successful 
the missile was in killing the target or the total probability 
of kill given a single shot (Fkss). 

An analysis of the final outcome "target killed" is useful 
for determining the sequential steps for the algorithm of 
engagement. 

The analysis of the final outcome "target not killed" is 
useful in pointing out the events that might cause a failure 
in killing the' target. The elements and events that lead to 
these two final outcomes should be considered if we want our 
model to be more realistic. In the real world, these elements 
and events can be countless. In our model they are limited to 
a few obvious ones, easy to understand and - varied enough to 
cover different areas. 
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A target not killed means that it was able to fly through 
all three envelopes successfully. An EEA should be conducted 
for each envelope, but since they differ only slightly from 
each other, only one EEA will be done for the long range 
envelope. The differences from the other two envelopes will be 
mentioned at the end of this chapter. 

The target may not be killed for the two following 
reasons : 

- the ship fires and misses 

- the ship does not fire at all 

If the ship fires and misses it will be for the following 
reasons : 

- missile malfunction (versus a failure to function) 

given a good fire control solution 

- fire control system malfunctions 

One category of missile malfunction could be a hardware 
problem, leading to a technical failure of the missile's 

guidance unit or the propulsion unit. Another category could 
be the guidance algorithms or the software was not 

sophisticated enough to prevent the missile from being 
diverted away from the target or simply an unsuitable 
navigation law. 

The fire control system malfunction can lead to a bad fire 
solution. For example, inaccurate parameters or wrong 
information about the target can be the cause. If there is a 
tracker problem or a successful use of the ECM by the enemy, 
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the target will not be tracked properly and a good solution 
will then be difficult to get. Another malfunction of the 
fire control system can be caused by a bad human— machine 
interface; the operator can be too quick or too slow to react, 
or he can push the wrong button, or becomes confused and 
misinterprets information. The consequence can be an 
inappropriate choice of the missile, such as a too slow 
missile for a fast or highly maneuverable target, or an 
inappropriate choice of the time of launch that may lead to a 
target out of range after firing the missile. Finally, a 
mechanical failure of a switch or a light to indicate a state 
(not ready or target not acquired, etc. ) can be the source of 
a malfunction. 

If the ship does not fire it could be for the following 
reasons. A fire control solution is generated, but the 
launcher fails to function or is not in a proper state or 
position. Also, a Command and Control problem may prevent the 
firing despite the availability of a solution. A safety 
factor or a management problem could be the source. 

If a fire control solution is not generated, the ship 
cannot fire. This can be caused by several different events: 

- target not detected or detected too late because of an 
ESM malfunction or sensors malfunction or inadequate 
perf ormances 

- target detected but not identified/designated properly 
because of the IFF malfunction or failure 

- target designated correctly but not assigned/acquired 
by the fire control system, and that can be caused by a 
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communication problem between the sensors and the FCS . 
A defective manual assign/acquire phase can be caused 
by the delay of the operator 

- target acquired but not tracked because of a tracker 
failure or an evasive action taken by the target (lost 
track ) etc . . , 

Figure II-l shows the relationship of these events. 

The main difference between the long range envelope and 
the medium range envelope is the detection-designation phase. 
Because that was accomplished in the previous envelope and 
needs to be done only once, there is no need for its 
repetition. The same applies for the short range envelope. 
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TARGET SURVIVES 
AND APPROACHES SHIP 
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Figure 1 1 -1 . Fault tree of the end-event "Target survives 
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Figure II-l. Fault tree of the end-event "Target survives". (Continued) 
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Figure II— 1 - Fault tree of the end-event "Target survives". (Continued) 



III. DESCRIPTION OF I SEAS 



A. BRIEF DESCRIPTION 

ISEAS is an interactive simulation of an engagement 
between a ship and an attacking aircraft. It simulates the 
functions carried on in the Combat Information Center, from 
the search phase to the eventual kill of the target. Three 
modules have been written separately by NFS students at the 
time of this thesis. These three modules have been combined 
by the author into one large program. These modules are the 
radar module dealing with the detection and the tracking 
phases; the weapons module, simulating the flyout of a missile 
from the launch to the intercept point; and finally the Combat 
Information Center module which solves the firepower equation, 
computes the geometric parameters of the engagement, and 
processes the information about the target and displays it in 
text and graphics modes. 

Each module is written in the C language and uses two 
types of variables. Local variables are used exclusively by 
their module where they are declared. Global variables are 
exchanged between modules so they need to be declared outside 
of the modules. In, this way, they are available for every 
module that needs them. Typical global variables are the 
target's position X, Y and 2. 
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These three modules are explained briefly below, pointing 
out their outputs and how they work with each other. 

B. RADAR MODULE 

The radar module [Ref. 1] is the first one called by the 
program when the simulation starts. Logically, it is the first 
phase of an engagement; the target has to be detected if any 
engagement is to occur. This module simulates the major 
functions of a radar, from the generation of an illuminating 
beam to the energy losses. to the detection of an obiect. The 
weather, any jamming by the enemy, and the characteristics ahd 
limitations of the radar are taken into account. This module 
simulates two different types of radars; the search radar and 
the tracking radar. The input for this module is the actual 
position of the target, which at this time is simply on a 
straight and level flight path. The output is the position of 
the target as tracked by a radar: contaminated by the 
uncerta ini ties and approximations because of the present 
techno 1 ogy . 

C. WEAPONS MODULE 

This module [Ref. 2] becomes active when the firing of a 
missile is ordered from the CIC. After the order, it reads the 
present position of the target and its predicted intercept 
point solved for by the CIC module. Then it calls the launcher 
function where the initial azimuth and elevation of the 
launcher are computed. Following is the call for the flyout 
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function. In this function the target's position is read 
every time it iterates, then the atmospheric conditions 
(density and temperature) are computed. After that the lift 
slope curve coefficient is calculated in the preparation for 
the solving of the angle of attack. The missile's thrust and 
weight are computed next. Finally, the guidance function is 
called and the angle of attack is solved. The output of the 
flyout function is the updated missile’s velocity and 
position. The updated position of the missile is sent to the 
CIC module. where it is processed and displayed. 

D. COMBAT INFORMATION CENTER MODULE 

This module is the controller of the program. Basically, 
it performs three different tasks. The first one is to solve 
the firepower equation, producing the firing geometry and 
analyzing it. The second one is to "orchestre" the different 
modules and subroutines of ISEAS on a time step event basis. 
The third one is to gather the information available and 
display it in a useful way. For details of how these tasks 
are carried out refer to Chapter IV in this thesis. 

E. FLOW OF THE PROGRAM 

Figure III-l shows the flow of the program. First, the 
input program is run and the variables are either changed or 
the default values are kept. These variables are saved in disk 
files. Then the simulation is run. It asks for which data 
files to be loaded. Then, the simulation starts by displaying 
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Figure III— 1. 



Flow of the Program 
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the tactical screen with the information on own ship only. 
The radar module is called and the search radar is activated. 
The search for an object starts into a loop. The detection of 
the target turns on a flag whithin the radar module when it 
reaches the maximum detection range. The later is computed 
after the technical characteristics of the radar and the 
weather conditions are checked. Once the detection flag is 
turned on, the target is reported to the CIC module where it 
is identified and designated as hostile. The firing geometry 
solved by the firepower equation and the closest point of 
approach (CPA) module is drawn on the tactical screen, Thus, 
up to now, we have a tactical screen where the maximum 
detection range, and the long, medium and short ranges of the 
missiles are drawn. On the 'text' screen, information 
concerning the radar being used and the weapons available are 
displayed. Also, information about the newly detected target 
is available. The CPA and its location are computed and drawn 
on the tactical screen along with the predicted flight path of 
the target. The position of the target is updated by 
continuously calling the radar module. This position is 
displayed on the 'text' screen and plotted on the tactical 
screen as it becomes available. Now, the user is looking at a 
radar-type screen with the necessary information permitting 
him to make a decision as to when to fire and which weapon to 
use. The CIC module, by this time, has already computed the 
possible firing sequence of the missiles depending on the 
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range of the target from the ship. The individual probability 
of trill given a single shot (Fkss) for each firing, and the 
total probability of kill are also available. Thus, the user 
knows what are his chances to kill that target. The plotting 
of the target’s position continues as long as it is within the 
detection range and the user doesn't take action. When he 
decides to fire at the target, the tracking radar is turned on 
by calling the tracking module. The later takes over the 
position updating with better accuracy. The user assigns the 
target to a weapon (long, medium and short range missiles 
provided they are available) and when he decides to fire, the 
weapons module is called, and the launcher position is 
computed and a ready signal is sent back to the CIC module. 
Now, the CIC module will be calling the tracking function, 
within the radar module, getting the instantaneous position of 
the target, displaying it and passing it to the weapons 
module. The latter goes to its flyout function and "flies” the 
missile accordingly and returns the missile’ sp os iti on to the 
CIC, which plots it on the same screen. This sequence is 
followed until the missile stops flying, then an analysis of 
the firing is done. The result is displayed in the form of a 
Pkss number and a total probability of kill if previous 
missiles were fired at that target. If the user decides to 
fire another missile, a delay time is observed, and the 
weapons module is called again in the same manner as before. 
The user keeps firing as he wishes provided some weapon is 
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available and ready to be fired, and as long as the target is 
within the missile's envelope. Once the target is killed 
(which is determined by a kill threshold value), or when it 
leaves all envelopes, a summary of the engagement is displayed 
as to how many missiles were fired and when, their intercept 
range and their Pkss. A total probability of kill is also 
displayed , The CIC module calls the radar module, setting the 
flags to their initial values, and activating the search radar 
for a new target. 

This simulation goes on until the user decides to stop it 
or until he expends all his missiles. 
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IV. 



THE DECISION MAKING 



A. DISCUSSION 

This chapter deals with the operational side; the 
decision-making and the analysis of the user's action. To 
define the steps needed for this task, a tactical situation is 
developed pointing out the events happening during an anti-air 
operation . 

TACTICAL SITUATION: The tactical situation may be best 

developed by listing the sequence of events which lead to the 
interception of the target by the surface to air missile. The 
following events would normally be expected to occur in 
sequence : 

1. DETECTION: The attacking' ai rcraf t is detected in some 
manner. The detection will be conducted from the ship 
with the surveillance radar. 

2. IDENTIFICATION: The target must be identified as 

friend or foe. If the target is identified as enemy, 
an alert must be given. 

3. EVALUATION: Information as to position and altitude, 
direction of motion, and probable type of attack is 
needed in order to direct the specific means of inter- 
ception that may be best employed at least risk to 

the ship. 

4. ASSIGNMENT: After the target has been evaluated, the 
the user will assign a weapon to it with the assistance 
of the command and control unit and the fire control 
system. 
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After receiving the information needed from the search 
radar and the tracking radar, and also from the weapons unit, 
the Command and Control unit will have to determine an 
optimised weapon assignment. The Fire Control system will 
generate a solution suitable to the target, and the Command 
and Control unit, given the tactical and management 
requirements, will decide on the engagement parameters, so the 
fire order will be issued at the right moment with the 
appropriate conditions. 

As is seen in the previous paragraph, two different units 
are involved so. far; the Command and Control unit, and the 
Fire Control system. 



B. COMMAND AND CONTROL UNIT 

Every element in the Combat Information Center, from the 
sensors to the weapons, performs differently in the various 
situations presented. For example a missile may have a better 
probability of hit when the target is closer. A sensor 
performs better when it is not functioning in a hostile 
environment. A target is easier to counter when the 
appropriate weapon is readily available and so on. The level 
of performance can be measured as a percentage of a maximum 
capability. Another way to think of it is how much can you 
expect that element to perform correctly. Usually, this 
performance is related to one or more parameters, such as a 
range, velocity, availability etc.. The value, which is a 
probability number, is either fixed or variable, found in a 
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table, guessed, or computed from a function. In this 
situation, a decision is more difficult to make because the 
decision maker is not judging if that element will function or 
not. Rather he will have to evaluate its performance with a 
value ranging from zero to one. Furthermore, a complete 
action involves many different elements. Consequently, the 
overall evaluation of the mission's result is more difficult 
to determine. The Command and Control unit plays a vital role 
in this process. With the appropriate set of rules foi 
computing pt obabi 1 i ties and finding expectations and 
sophisticated tactics and algorithms, the Command and Control 
unit can provide the decision maker with an optimized reaction 
in any tactical situation. The advantage of the Command and 
Control unit, which is merely a processor loaded with a 
sophisticated program, is speed and memory. It is built to 
deal with various tactical situations, to interact with the 
different sensors and equipment on board, to process 
information according to a set of rules at very high speed, 
and to recommend an optimal action in a easy way to use. 
Usually, the output is expressed by means of graphics on 
tactical CRT screens. 

The Command and Control unit implemented in ISEAS gives a 
simple example of how a number of missiles will be suggested 
to fire at given times so they will achieve a wanted total 
kill probability. Also it gives the individual and the 
accumulated kill probabilities for each missile fired or to be 
f i red . 
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Another use of the Command and Control unit is to help 
manage the missiles regarding the availability and the kind of 
mission and the type of target to deal .with. There will be 
missiles on board for different ranges (long, medium and 
short) with different capabilities; mainly speed and 
navigation law. Once the information about the target is 
available, the Command and Control unit will help in choosing 
the right weapon for the target. The problem of managing the 
number of missiles could be solved by the same unit. The 
operator can either fire at the target every time it is 
possible regardless of how many missiles will be expanded, or 
wait for the target to get closer so the expected kill 
probability reaches a preestablished value and then fire. In 
the later case less missiles would be fired and the Command 
and Control unit is appropriate for such computations. 

The Command and Control unit module in ISEAS can be 
modified and improved to deal with more complex problems and 
take into account more parameters. 

The first problem that this unit solves in this program is 
the performance of the missile measured against the intercept 
range. A missile is supposed to have a maximum range and a 
minimum range. In this model, the missile will have a high 
probability of kill at the minimum range ( RN ) and a low 
probability of kill at the maximum range ( RX ) and a zero 
probability outside its operating range (either beyond the 
maximum range or within the minimum range ). The probability 
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of kill given a single shot (Pkss) will be a function of the 
intercept range, varying between two nominal values within the 
operating range. The simplest type of function that will do 
this task is a linear function (Figure IV-1 ) . The function 
will have the following form: 

Pkss (range) = A * range + B ; RN < range < RX 
Pkss( range) = 0 if' range > RX or range < RN 
This .function can be used in two different ways; to 
compute the kill probability given an intercept range, or find 
the range at which a specific probability of kill can be 
achieved. 



Pkss 




RANGE 



Figure IV-1. Linear function of Pkss vs. Range 
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This function gives the probability of kill given a single 
shot, which is considered a success p versus a failure q. The 
accumulated probability of kill for multiple shots is 
given in the following form: 

Pk(total) =1 - q ** n 

with q = 1 - Pkss for each missile fired. 

But since the value of Pkss varies for each missile 

(because each missile is fired at a different distance) the 

formula for the total probability of kill will be 

Pk(total) = 1 - (1 - Pkss(i)) where Pkss(i) is 

the probability of target kill for the ith missile fired. 

Another use of the function of the probability of kill is 
to determine the intercept range at which a given value of the 
kill probability is sought. 

This is done by inverting the initial function to the 

foil owing 

range = (Pk(range) - B) / A . 

This range is the intercept range. For the intercept to 
occur at this range the fire control system has to solve the 
geometry and find the firing range which is prior to the 
intercept range. 

C. FIRE CONTROL SYSTEM 

One of the main functions of the fire control system is to 
determine how many missiles can be fired at an incoming 

threat. The solution is given by the FIREPOWER equation 
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derived here. The raid could be either an aircraft or a 
missile, with the following conditions: 

- the raid is moving at a constant velocity, constant 
direction and constant altitude (non maneuvering-target ) 

- the heading is not necessarily toward the ship 

- the geometry and the equations are treated in three 
dimensi ons 

- the missile has a maximum intercept range and a minimum 
firing range 

To fire the maximum number of missiles, the first missile 
must be fired so that it will intercept exactly at its maximum 

range. Subsequently missiles must be fired as fast as 
possible, until the raid has closed within the minimum firing 
range or left the envelope. In developing this equation, a 
shoot-look-shoot policy has been used with a variable look 
delay time. 

The shoot-1 ook— shoot model is especially valid when only 
one fire control radar or guidance system is available for 
that particular target. To solve the firepower equation, the 
closest point of approach (CPA) is needed. A module called CPA 
solves this problem. ■ Once the CPA is solved for, the geometry 
for the firepower equation can be drawn and solved. 

1 . The C 1 osest Point of Approach Modu 1 e 
The variables used in this module are: 

CPA : the closest point of approach in the X-Y plan. 

TCPA : the closest point of approach in the X-Y-Z 

space . 

XT.YT.ZT : the initial position of the target. 
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VT , CT 
XT , YT 
XS , YS 
VS ,CS 



: the target's velocity and course. 
: the target's velocity components. 
: the initial position of the ship. 
: the ship's velocity and course. 



XS.YS : the ship’s velocity components. 

X,Y : the components of the target's relative 

velocity. 

The output will be the true TCPA (three dimensions). 
The geometry is shown in Figure IV— 2. Zero degree heading is 

considered along the X— axis and positive to the left 

(Cartesian coordinates). 

The ship to target horizontal distance is computed 
knowing the initial positions. 

R = SQR [(XS-XT)**2 + (YS-YT)**2] 

The angle A from the target to the ship in the horizontal plan 
is given by 



A = ARCTAN [(YT-YS) / ( XT-XS ) ] 

'The components of the target’s relative velocity are 
computed next. The ship is taken as the reference so that 

X = XT - XS 
Y = tfT - YS 

where 

XS = VS * cos(CS) 

YS = VS * sin(CS) 

XT = VT * cos(CT) 

YT = VT * sin(CT) 
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>- 



Figure IV-2. 



Closest Point of 



Approach 
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(SA'SX) dIHS 



The angle of the relative velocity is found from 
B = ARCTAN (Y/X) 

A new angle C is computed from the geometric 

difference of the two previous angles. This angle gives the 
geometry of the final triangle that will lead to the CPA in 
the X-Y plan. 

C = A - B 

The CPA in the X-Y plan is the opposite side of the angle C. 
Thus . 

CPA = R * sin ( C ) 

The true CPA takes the altitude of the target into account. 

TCPA = SQR [ ( CPA ) **2 + (ZT)**2] 

TCPA will be referred to as the CPA meaning the true CPA. 

A test is made on the angle C. If it is larger than 90 
degrees (absolute value) no CPA case exists. Geometrically 
speaking, the target will not get any closer than where it is 
at the time of the computation. 

2 . The Fi repower Equat i on Modu 1 e 

Now that we have the CPA we can construct the geometry 
for the engagement. For the purpose of clear demonstration two 
different geometries will be considered: 

CASE I: the target is on any heading (variable CPA), 

but with no delay look time (T = 0). [Ref. 3] 

CASE II: this is the general case, the target is on 

any heading, with a variable delay look time T. 
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The following terms, in addition to the ones already 


defined 


in the CPA module will be used in the derivation: 


SAM : 


surface to air missile. 


RX : 


maximum intercept range of SAM. 


RN : 


minimum intercept range of SAM. 


VM : 


speed of SAM. 


t : • 


time. 


n : 


number of missiles after initial intercept at RX . 


TA- : 


the target angle, from the target, between the 

target's flight path and the line of sight. 


LA : 


lead angle, from the ship .between the missile's 
flight path and the line of sight. 


RT : 


the distance traveled by the target between 
successives intercepts. 


RI 


the intercept range, measured from the ship. 

As it was stated before, the first missile must be 



fired so that it will intercept at the maximum range of the 
missile. The equation will be derived for the additional 
number of missiles fired after the raid has reached P.X , and 
one missile will be added to this number. 

CASE I: target with any heading and no delay time. 

An equation giving the number of missiles that can be 
shot at an incoming threat as a function of the target's 
closest point of approach is derived below. The geometry 
involved in this model is shown in Figure IV-3. 
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From this figure it can be seen that the smaller the 
CPA is, the more missiles that can be fired. The flight path 
of the target within the missile's envelope increases as the 
CPA decreases. 

The number of missiles that the ship can fire during 
the time the target is within range is a function of the 
parameters RX , RN , VM , VT and the missile's flight geometry. 
Figure IV-4 illustrates how the flight geometry changes with 
two different laws of navigation (pursuit and lead angle), 
thus changing the time of flight and the distance to the 
intercept point. 

In the following derivation, the lead angle navigation 
with a constant line of sight is chosen for three reasons: 1) 
this course is more efficient in terms of the SAM time of 
flight before it reaches the target; 2) the equation is easier 
to develop; 3) it is often used by many missiles. 

Figure IV-5 illustrates the problem to be solved in 
calculating the number of possible shots as a function of the 
CPA. The initial target angle TA can be defined as 

TA = Arcsin ( CPA/RX ) 

The lead angle is derived using Figure IV-6. For an 
intercept course to exist, the components perpendicular to the 
line of sight of both the target speed (Vta) and the missile 
speed (Vma) must be equal. Thus 

Vta = V T * sinTA 
Vma = VM * sinLA. 
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LEAD ANGLE NAVIGATION 




Figure IV-4. Flight geometry changes with two different 
laws of navigation. 



39 



II 



I— 

u 



c 




40 



Figure IV-5. Geometry of the engagement with no delay time (Case I). 



TARGET 




Figure IV-6. Lead angle and target angle. 



therefore 

VT * sinTA * VM * sinLA 

and hence 

LA = Arcsin [(VT * sinTA) / VM ] 

Once the lead angle equation is defined, the problem 



of Figure 


IV— 5 an be 


solved 


for 


RT and 


RI . Knowing 


both 


lead 


angle and 


target 


angle , 


the 


third 


angle ( PHI 


) 


of 


the 


triangle 


described 


by RX, 


RT 


and RI 


in Figure IV- 


-5 


can 


be 



determined by using 

PHI = PI - (TA + LA) 

RT and RI can be found by the law of sines 

sinPHI / RX = sinTA / RI = sinLA / RT 
RI = RX sinTA / sinPHI 
RT = RX sinLA / sinPHI 
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If a shoot-look-shoot firing doctrine is used, and a 
zero delay time between shots is assumed, then the intercept 
range for the previous shot becomes the firing range for the 
succeeding shot, and the target angle for the succeeding shot 
can be expressed as 

TA(n) = Arcsin[CPA / Ri(n-l)] 
for a closing target, and 

TA(n) = PI - Arcsin [Rcpa/Ri ( n— 1 ) ] 
for an opening target. 

Another way to compute this angle regardless of the 
target's situation is to use the angles of the previous firing 

TA ( n ) = PHI(n-l) + LA(n-l) 

Now the initial lead angle LA can be calculated, and 
RT(2) and RI(2) calculated, and so on, until the target has 
flown out of SAM range. 

To get the maximum number of missiles, a missile must 
be fired to intercept exactly at the maximum SAM range. 
Therefore, each intercept shown in Figure IV— 3 is in addition 
to the initial intercept at maximum range. 

The computations of the firing parameters, the 
intercept range mainly, are iterative. Each firing geometry 
depends on the previous firing, and so on until the initial 
intercept which is RX . As long as RI is smaller than RX the 
iteration will be executed, and the number of missiles that 
can be fired is incremented. When this condition is not true 
anymore, the iteration is abondonned and the total number of 



42 



missiles to be fired is n plus the first missile for the 
initial interception. 

CASE II: target with any heading and a delay time 
For this general case, the target is approaching the 
ship from any direction on any course as in case 1, However, a 
delay time T is adopted for observation of the missile at the 
intercept point to see if another missile is needed or not to 
kill the target. In the previous case, the intercept point was 
considered to be the initial firing range for the succeeding 
missile. In this case, the target will travel a certain 
distance corresponding to the delay time T before a new 
engagement geometry will be solved for. This distance is the 
same for all the shots since the delay time and the target's 
velocity are fixed. To solve the firepower equation the target 
angle (TA) after the delay time is required, are needed first. 
Three new parameters are defined in Figure IV— 7 : 

- NFHI is the angle between the flight path of the target 
and the previous intercept range RI . 

- NL is the angle opposite the distance DIST corresponding 
to the delay time. 




Figure IV— 7. Geometry of the engagement with delayed 
time . 
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- NR is the new initial firing range used to compute the 
target angle. 

The equation for TA is solved in the same manner as 
before, but there are more steps in this case. They are: 

NPHI = Arcsin [CPA / RI ] 

BIST = T * VT 

NR = SQR { DIST ** 2 + RI ** 2 - 2*DIST*RI*CosNPHn 
NL = Arcsin [DIST * CPA / (NR * RI ) ] 

TA = NPHI + NL 

With TA determined, the computation is continued as 
for the previous case to solve for the RI . 

This iteration tests for the next intercept range to 
see if it is valid to fire another missile or not by comparing 
RI to RX. 

When the intercept of the target is possible but the 
target itself is far from the ship, it takes a relatively long 
time for the missile and the target to reach the intercept 
point. Thus, very few missiles will be expended. But the 
closer the target is from the ship the shorter is the time 
needed to reach the intercept point. Thus, many missiles will 
be expended within a few seconds. 

Since many types of missiles have a minimum intercept 
range, the iteration solving the firepower equation is 
improved to take the minimum range RN into account. 
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When the target closes within the minimum range the 



iteration stops and then restarts immediately 
intercept range at RN for the departing target, 
then continues for the departing target as before 



with a new 
The iteration 
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V. 



THE INPUT MODULE 



This chapter describes the opening screens of ISEAS and 
the input screens which give the user the choice to change any 
of the variables, 

A. DIFFERENT SCREENS 

These screens are all in graphics mode (versus system's 
text). The graphics library used in ISEAS is HALO, which is 
considered to be the standard for the Graphics Kernel System 
( GKS ) . For detailed information about the use of HALO in 
ISEAS see Chapter VI. (Graphics). 

1 . The Opening Screens 

When ISEAS is executed, the first screens displayed 
are the opening screen (Figure V— 1 ) and the explanatory 
screen. Their function is to inform the user of ISEAS, what 
does the word ISEAS stands for and its developers. The 
explanatory screen briefly describes ISEAS, discussing what it 
does, what is its scope, how to use it, and where to refer to 
for more information. 

2 . The Input Screens 

As it was described in Chapter III. , ISEAS is a program 
that consists of three different modules; the radar module, 
the weapons module, and the combat information center module. 
For versatile use of ISEAS, the program has an input 
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Figure V-l. The Opening Screen, 
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capability in which the user can change some or all of the 
variables and judge their influence on the outcome of the 
simulation . 

In ISEAS, every module has an one or more independent 
input screens such as shown in Figure V-2 . Basically, they 
are all built the same way. First, since they are graphics 
screen', HALO’S environment is initialised and a viewport (i.e 
window) covering 75% of the screen is set(the other 25% is 
reserved for the scanning of the values of the variables as 
exp 1 ained later). This viewport is filled with a color as a 
background color and mapped from 0 to 100 on the horizontal, 
and 0 to 50 on the vertical with the origin (0,0) located at 
lower left corner of the display monitor. This mapping is 
necessary so the application can be device independent, and 
also for the locating of HALO's text cursor. On this viewport, 
a title for that screen is displayed. Since this is the input 
module, the title for all the screens is "INPUT MODULE". A 
subtitle displayed below it informs the user as to which set 
of variables are actually being ‘ input. There are four 
subtitles and thus four input screens; the search radar, the 
tracking radar, the weapons, and the combat information center 
input screens. On each input screen the list of variable names 
that can be changed is displayed. Each variable name bears a 
letter (from a to z) as an individual identifier. Note that 
the values of these variables are not displayed yet. They are 



48 



Figure 




- 2 . 



Search Radar input screen with no values 
display yet. 
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stored in disk files which are identified by filenames that 



match their module name; the extension of these files is "DFL" 
for default. The program scans the variables values from these 
files and automatically displays them next to their variables 
names respectively (Figure V— 3 ) , The user is then asked if he 
wants to change some of the variables. A message at the top of 
the screen refers the user to the manual (to be prepared) for 
the meaning and limits of these variables. If the user chooses 
to change some variables and presses the character 'Y' , the 
program switches to the "input mode", This is a 'while loop' 
with the condition being i (the input character) equal to 'Y'. 
A SWITCH instruction branches to the appropriate CASE 
statement based on the letter entered and identifies the 
variable name to be changed. The new value is entered at the 
keyboard. The old value is erased then (Figure V-4 ) , and the 
value newly entered is displayed. As long as i is equal to 'Y' 
the user is in the input mode which allows him to change any 
variable as many times as he whiches and in any order he 
wants. He can leave the input mode by pressing the character 
' x' which results in the branching to the CASE * x' where i is 
set equal to ’N’. The condition for the while loop is not true 
anymore and the program exits the input mode. The display 
of the values on these input screens was not easy to create. 
HALO has no capabilities of inputing or outputing floating 
numbers from the console or the keyboard and no text can be 
sent directly from the keyboard. It displays text (not numbers 
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Figure V-3 . 



Search Radar input screen with the default 
values displayed. 
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Figure V— 4. 



Input screen with a value [s.BW] to be 
changed. 
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such intergers or float values) that was stored in a string. 
Thus all the variables values are converted from integer or 
float numbers to character text! This was done with the help 
of the Lattice C compiler, except that scanned variables are 
echoed on the display monitor, which destroys the graphics on 
the viewport! The method used to fix this inconvenience was 
the allocation of 25% of the total area of the display for 
text input and number input. The rest was allocated to the 
graphics with the viewport defining the limits between the two 
modes . 

When the user exits the input mode he is asked if he 
wishes to save his set of data in a disk file (for future use, 
reference or comparison) (Figure V— 5 ) . In the affirmative, he 
is asked to enter the filename and extension of his file. 
Since it is a disk file the filename can be no longer than 
eight characters, and the extension no longer than three 
characters. After that, the program opens a disk file by that 
name and stores the data in it. This last action ends the 
current input screen, and the user is prompted to hit -any key 
to continue. A new input screen is displayed with the 
corresponding title and variables list, and so on, until all 
the input screens are displayed. 

ISEAS consists of two separate programs. The first 
program, the input program, is used to change variables and 
to create data files. The second program is the actual 
simulation. The later asks for a data file to load. The user 
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Figure V-5. An input screen with the option to save 
the data. 
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will either choose the default files or his own, which was 
created by the first program. The default files are provided 
with the program. But if they are not available for a reason 
of a loss or other reason, they can be generated easily. The 
variables will be zeroed out because of their nonavailability 
(Figure V— 6 ) but they can be changed to new default values and 
stored in the default files. The same procedure can be 
followed if the default files need to be changed. 

B. WHAT ARE THE OPTIONS? 

The set of variables used in ISEAS is divided into three 
parts according to the three different modules. In the 
following, these parts are described as to how many screens 
each part is using and the list of variables displayed on that 
screen with the definition of each variable. 

1. RADAR INPUT 

This module has two input screens; the search radar 
screen and the tracking radar screen. 

a. Search radar screen (Figure V— 3 ) 

This screen has twenty variables, their default 
values are stored in the file called "SEARCH . DFL “ . Their 
definitions, units and default values are: 

a . P J : jammer power, 0.1 Watts. 

b. RT : recovery time, 0 . 5E-6 seconds. 

c. BAZ : azimuth beam width, 1.5 degrees. 

d. GJJ : jammer antenna gain, 0.0. 

e. T : pulse width, 2E-6 seconds. 
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Figure V-6. An input screen with the variables 
"zeroed out" 
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f . BEL 



f . BEL : 


elevation beam width, 10.0 degrees. 


g . JR : 


jammer to radar distance, 45.0 meters. 


h . AZR : 


antenna scan rate, 120,0. 


i . POWER : 


transmitter power, 3 . 0E5 Watts. 


j . FAN : 


false alarm number, 1 , 0E10 ? 


k . NF : 


noise figure, 5.0 ? 


1 , PL : 


plumbing loss, 2.0 ( dB ?) 


m . NS 


number of scans, 3. 


n . BL 


beam loss, 3.0 (dB ?). 


O . EA : 


effective area, 0.8 square meters. 


p . AA : 


antenna altitude, 30 meters. 


q. ST : 


system temperature, 290.0 degrees Kelvin. 


r . PRF : 


pulse repetition frequency, 1000.0 Herts. 


S. BW : 


receiver bandwith, 5.0E6 Hertz. 


t. FF : 


transmitting frequency, 9.0E9 Hertz. 


X.EXIT : 


quits the input mode. 



b. Tracking radar screen (Figure V— 7 ) 

This screen has twelve variables. The variables 
names are similar to those of the search radar. However, to 
differentiate them from the previous ones they all have the 
letter T at the end (for Tracker). The values are different, 
and some variables found in the search radar screen don’t 
exist in the tracking radar screen. These variables are stored 
in a file named "TRACK. DFL". For more details refer to Ref. 1. 
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Figure V-7. 



Tracking Radar 



Input Screen. 
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CIC INPUT 



2 . 

This module has one input screen. Its variables are of 
the ship's situation, and the CIC. The variable values are 
stored in the file called "CIC.DFL". 

a . CS : ship’s course, 0.0 degrees. 

b. VS: ship's speed, 25,0 Km/h . 

c. T : time increment for the simulation, 0.1 seconds. 

d. LT: look time for the look-shoot-look policy, 4 

seconds. 

e. TK: treshold for deciding a kill has occured, 99.5% . 

x. EXIT. 

3. WEAPONS INPUT 

This module was not available at the time of writing 
this thesis. 
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VI . GRAPHICS 



As was discussed in the Introduction, the graphics 
represents a large part of this thesis. One of the best ways 
to take full advantage of the stream of information flowing in 
the CIC is to display it in the form of graphics on tactical 
display CRTs. 

In ISEAS, the critical information is gathered from other 
modules, processed and displayed graphically on the screen, 
simulating what is actually happening or what might happen in 
the near future. Thus, the user is saved a big step of 
imagining in his mind, the scenario of the engagement. 

A. HALO 

The graphics portion of ISEAS is written in the C 
language with the graphics library HALO. The later follows the 
Graphics Kernel System ( GKS ) . HALO, version 2.30, and lattice 
C, version 3.00G, are the actual library and compiler used. 

HALO is a library of subroutines offering the programmer 
the possibility of adding graphics to his application program. 
In the following paragraphs, some concepts used in ISEAS are 
explained permitting the user to understand how the simulation 
is conducted and also permitting an improvement of ISEAS. 

HA.LO has first to be initiated.. During this initiation 
the interfacing devices used have to be defined, especially 
the type of monitor used. Once the graphics mode is initiated 
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and the appropriate device driver is loaded, graphics are then 
possible to draw. The actual version of ISEAS is designed to 
run with either the IBM Color Graphics Adaptor (CGA) or the 
IBM Enhanced Graphics Aaptor (EGA). Any modification done so 
that ISEAS can run on another type of monitor can be easily 
done in the source code. 

ISEAS is fairly device independent. To achieve this, the 
world coordinates are used throughout the program. The world 
coordinates (WC) are a set of coordinates, user defined, that 
map the screen to the range of coordinates chosen by the 
programmer, A point is located with these coordinates 

regardless of the resolution of the monitor used. In ISEAS, 
the world coordinates range from zero to one hundred on each 
axis for all of the graphics screens. The viewport (or window) 

is the part of the monitor dedicated to graphics. The whole 

surface of the monitor could be dedicated to one graphics 
screen, the limits of the viewport are then (0. 0,1.0) on each 
axis (i.e. , the whole length of each side is used as limit for 
the viewport). Some screens in ISEAS use the whole -display 
screen for one graphics screen. But mostly two viewports share 
the same display screen. The reasons for this partition are: 

a. two different viewports are needed, one for the 

graphics of the simulation and the other for the display 
of the information in text form. 

b. the surface of the monitor is not square. Thus, 
anytime a location is addressed, its coordinates have 
to be scaled to the aspect ratio. Otherwise, the 
range on the horizontal axis will be larger than the 
vertical axis. This will not fit the circular view given 
by the radar. 
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Because of this, a square viewport used by the graphics 
is located to the left and it ranges from 0.0 to 0,6325 on the 
horizontal axis (x-axis) and 0.0 to 1.0 on the vertical axis 
(y-axis), (0.0,0. 0) being the lower left corner. 

Non-destructive graphics are wanted most of the time. 
This is when the background or whatever was previously 
displayed is restored once the present image is 'erased' . 
Erasing an image is merely rewriting it twice. The second time 
it is written, the image is XORED with the first one and the 
background is restored. This is achieved by turning the XOF. 
mode on. Animation is also realized by the same procedure. A 
figure is drawn once, so that it appears on the screen, then 
it is drawn a second time so that it is 'erased' , then the 
coordinates of the figure are shifted toward the direction of 
the motion and redrawn and so on. One problem arises when 
using this mode with a colored background; the colors change 
and some of them don't show at all. The background used in 
ISEAS is light blue and the colors used are carefully chosen 
so that the user sees the screen (and the simulation and data) 
easily. 

B. THE SIMULATION SCREEN 

The graphics are used mainly to show what is physically 
happening during the engagement. Also, possible future 
engagement geometries, intercept points and range limitations 
are shown on the graphics screen. The preengagement phase and 
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the engagement itself can be followed on the different types 
of screens. The first one is a radar-type, display, where a 
circular view is adopted with the ship permanently located in 
the center. The second screen is a TV-type display (not 
discussed here), where the engagement is followed on the three 
different planes. Both of them have a viewport reserved for 
text, messages and data displaying. In the following these 
types of graphics screen are explained in details. 

The Radar-type display is the first screen displayed 
(Figure VI-1), and it is used for the general information 
about the tactical situation and the ship's situation. It has 
two viewports, a square one on the left that simulates a radar 
scope display, a second one to the right where messages and 
text information are displayed periodically. 

On the first viewport, the ship is located in the center. 
A pair of orthogonal axes are drawn with marks along them. 
Then four concentric circles are displayed. The outer one is 
fixed for all the simulations, and it represents the maximum 
detection range of the search, radar used, it is labeled ”DE". 
The other three circles represent the ranges for the long 
range, medium range and short range missiles on board. Mo 
other information is displayed before the detection of a 
target , 
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>*■!* 






Figure VI-1. Tactical screen with no contact 
made yet. 
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The second viewport (Figure VI-1) is divided into three 
windows showing- information about the radar used, the missiles 
and the target. The first window displays information about 
the radar; the search radar or tracking radar used, then the 
operating frequency, then the maximum detection range, and 
finally a message about the contact. At this stage, the 
message "NO CONTACT" is displayed. The second window contains 
information about the missiles. Their ranges, available number 
for each one and status (ready/not ready) are made known to 
the user. The third window, and the most important, gives 
informations about the target once it is detected. Before its 
detection, only the titles of the parameters are displayed 
with no data available. These parameters are the identity of 
the contact (ID) designating it as friend, enemy or unknown. 
Then the present position of the contact from the ship is 
displayed. The CPA and its position are also displayed. 
Finally, the speed and course of the contact are shown. The 
current target is non-maneuvering, thus the only information 
to be updated constantly in this window is the position and 
the range. When a contact is made, and that happens only when 
the target is within the detection range, its position and 
course are plotted instantaneously, then the CPA is drawn 
showing its location relative to the ship. On the 'text' 
viewport first window, the message "NO CONTACT" is replaced by 
a blinking highlighted message "WARNING! CONTACT". The third 
window displays all the information available about 
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the contact and updates its position while it is plotted on 
the 'graphics' viewport simultaneously. The target keeps 
flying, and the user takes action. The choices are discussed 
later, but if the target leaves the detection envelope, the 
updating of the position and the warning signal stop. A new 
message "CONTACT LOST!" is displayed, and the search continues 
for a new contact. 

During the simulation of the target's flight, the CIC 
module solves the firepower equation for the geometrically 
possible intercept points for the three different range 
missiles. Then these intercept points are displayed on the 
screen. The combination of the circles for the three ranges 
and these intercept points gives a complete idea to the user 
of what are his capabilities regarding that specific target. 
The user has the possibility to fire the missile of his choice 
at a time he judges appropriate. His judgement will be based 
on what he sees on the screen regarding the present tactical 
situation and his fire power capabilities. And that is the 
main idea behind using the graphics in the CIC. 
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VII. 



SUMMARY 



A computer program simulating an anti-air operation 
conducted from the C.I.C. of a ship was written in the C 
language to run on an MSDQS personal computer. This program 
simulates teh main functions of a C.I.C. and is incorporated 
into the NPS interactive simulatiron of an engagement at sea 
( ISEA.S ) with a radar module and a weapons module developed by 
others. The contribution of this thesis to ISEAS is weapons 
direction, decisi on— making and graphics display of the 
tactical information in a useful form. 

The present version of ISEAS is the product of a first 
iteration which was dedicated mostly to define the frame work 
and scope. It is intended to demonstrate basic features found 
in any CIC, and encounter between a ship and an attacking 
aircraf t . 

The author hopes that more students will take this task 
over and improve the present version and enhance it to a 
higher level with more capabilities. It is a task that can be 
carried to a high level of expertise. 
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APPENDIX A 



#include"stdio. h" 
#inc 1 ude "s tdl ib . h " 



float xl = 
float x2 = 
int black 
int blue 
int green 
int red 
int yellow 
int white 



0.0 , yl 
100. 0 ,y2 
= 0 ; 

= 3 ; 

= 2 ; 

= 12 ; 

= 14 ; 

= 15 ; 



0.0 ; 
50.0 ; 



int purple = 5 ; 
int liblue = 3 ; 

float PRF , AA, T. RT , EA , FF , PL, BL , BW , ST. NF , 

POWER, AZR . NS, BAZ . BEL, FAN, P J , G JJ , JR ; 
float AAT , EAT , BWT , CST , POWT , FFT , GLT , BLT , STT , NFT , T AZ , TBEL 



float cx , cy ; 
int h , w , p , m ; 
int i ; 
int pixels ; 
int switsh = 1 ; 

float vxl =0 . 25 , vy 1=0 . 0 , vx2 =1 . 0 , vy2=l . 0 ; 

char s [20 ] ; 

int border = -1 , back = -1 ; 



main ( ) 

int mode ; 

setdev ( "HALOIBME. DEV " ) ; 

mode = 4 ; 

initgraphics( &mode ) ; 

setwor 1 d ( &xl , &y 1 , Xx2 , &y2 ) ; 

Opening ( ) ; 

intro ( ) ; 

inputl ( ) ; 

changel ( ) ; 

input2 ( ) ; 

change2 ( ) ; 

c 1 oseg raphics ( ) ; 

1 
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Opening ( ) 

{ 

float cx , cy ; 
int h , w , p , m ; 
int i ; 
int pixels ; 

pixels =13 ; 

set 1 nwidth ( .Jpixe 1 s ) ; 

setcolor(Xred) 

box ( &xl , &yl , S.x 2 , &y2 ) ; 

cx = 50 . 0 ; 

cy = 49 . 0 ; 

movabs ( &cx , &cy ) ; 

f 1 ood( &1 ibl ue ) ; 

h = 2 ; 

w = 2 ; 

P = 0 ; 
m = 0 ; 

set text ( &h , &w , &p , &m) ; 

settextc 1 r ( &bl ue , Xbl acK ). ; 

cx = 40.0 ; 

cy = 44.5 ; 

movtcurabs ( &cx , Scy ) ; 

text( "WELCOME" ) ; 

cx = 40.0 ; 

cy = 41.0 ; 

movtcurabs ( &cx , &cy ) ; 

text ( "TO THE”) 

cx = 18.0 ; 

cy = 36.0 ; 

settextcl r ( &black , &black ) ; 

movtcurabs ( Xcx , &cy ) ; 

text( "NAVAL POSTGRADUATE SCHOOL") 
h = 3 ; 

settex t ( &h , & m , 2p , &m) ; 

settextcl r ( &red , ablack ) ; 

cx = 38.0 ; 

cy = 30.0 ; 

movtcu rabs ( &cx , &cy ) ; 

text ( "I") ; 

cy = 25 . 0 

movtcurabs! &cx , &cy ) ; 

text ( "S " ) ; 

cy = 20.0 ; 
movtcurabs ( Sex . &cy ) ; 

text ( "E" ) ; 
cy = 15.0 ; 
mov tcurabs ( Sex , Scy ) ; 

text ("A") ; 

cy = 10.0 ; 
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movtcurabs ( &cx , &cy ) ; 

text ( "S ’’ ) ; 

h = 2 ; 

settext ( &h , &w , &p , &m ) ; 

settex tc 1 r ( &ye 1 1 ow , &bl acK ) ; 

cx = 41.0 ; 

cy = 3 0 . 2 ; 

movtcurabs ( Sex , &cy ) ; 

text ( " INTERACTIVE" ) ; 

cy = 2 5.2 ; 

movtcurabs ( Sex , Scy ) ; 

text( "IMULATION of an") ; 

cy =20.2 ;• 

movtcurabs ( Sex , Scy ) ; 

text( "NGAGEMENT" ) ; 

cy = 1 5 . 2 ; 

movtcurabs ( sex , &cy ) ; 

text ( “T " ) ; 

cy = 10 . 2 ; 

movtcurabs ( Sex , Scy ) ; 

text ( "EA " ) ; 

h = 1 ; 
w = 1 ; 

settext ( Sh , &w , &p , &m ) ; 

settextc 1 r ( &bl ack , Sblack ) ; 

cx = 70 . 0 ; 

cy = 1.0 ; 

movtcurabs ( Sex , Scy ) ; 

text("Hit any Key to continue") 

i = getch ( ) ; 

setcol or ( sblack ) ; 

cl r ( ) ; 

return ; 

1 
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intr o ( ) 

{ 

float cx , cy ; 
int h , w , p , m ; 
inti ; 
int pixels ; 

pixels = 13 : 

set 1 nwidth ( &pixe 1 s ) ; 

setco 1 or ( Sred ) ; 

box ( Sxl , &yl , ?,x2 , $y2 ) : 

cx = 50. 0 ; 

cy = 49.0 ; 

mouabs ( &cx , Scy ) ; 

f 1 cod ( & 1 ibl ue ) ; 

h = 3 ; w = 2 ; 

p = 0 ; m = 0 ; 

settext ( &h , &w , &p , &m ) ; 

settextcl r ( &bl ue , Sblack ) ; 

cx = 34.0 ; cy = 42.0 ; 

movtcurabs ( Sex , &cy ) ; 

text ( •'INTRODUCTION" ) ; 

h = 1 ; w = 1 ; 

settext ( &h , aw , &p , am ) ; 

settex tc 1 r ( Spurpl e , £bl ack ) ; 

cx = 14,0 ; cy = 36.0 ; 

movtcurabs ( &cx , &cy ) ; 

text("This program simulates a Combat Information Center in 
action " ) ; 

cx = 10.0 ; cy = 34,0 ; 
movtcurabs ( &cx , &cy ) ; 

text("in an Anti-Air operation , A target flying at a 
constant velocity") ; 
cy = 32,0 ; 
movtcurabs ( Sex , Scy ) ; 

text("and at a fixed altitude is generated. The ship has for 
mission the") ; 
cy = 30 , 0 ; 
movtcu rabs ( acx , S,cy ) ; 

text ( "counter ing of this target and ultimately its 
destruction . " ) ; 

cy = 22.0 ; 
movtcurabs ( Sex , Scy ) ; 

text ( " THIS PAGE WILL BE RESERVED FOR A SHORT INTRODUCTION 
TO THE PROGRAM"); 

h = 1 ; 

W = 1 ; 

set text ( &h , , .?p , jm ) ; 

se t tex tc 1 r ( sbl ack , Sbl ack ) ; 

cx = 70.0 ; 
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cy = 1.0 ; 

movtcurabs ( &cx , &cy ) ; 

text("Hit any key -to continue") ; 

i = getch ( ) ; 

setcol or ( Sblack ) ; 

clr() ; 
return ; 

} 
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inputl ( ) 

I 

FILE *pointerl ; 

setviewport( avxl , avyl , avx2 , avy2 , aborder , aback ) ; 

pixels = 13 ; 

setl nwidth ( apixels ) ; 

setcol or ( ared ) ; 

box (Sxl.&yl, ax 2 , &y2 ) ; 
cx = 50.0 ; 

c y = 49.0 ; 

movabs ( acx , acy ) ; 

f 1 ood ( a 1 ibl ue ) ; 

setxor ( aswi tsh ) ; 

h = 3 ; w = 2 ; 
p = 0 ; m = 0 ; 
settex t ( &h , &w , &p , am ) ; 

set tex tc 1 r ( agreen , abl ack ) ; 

cx = 30.5 ; cy = 45.0 ; 
movtcurabs ( acx , &cy ) ; 

text ("INPUT MODULE") ; 
settextcl r ( &red , ablack ) ; 

h = 2 ; m = 1 ; 

settext ( &h , &w , &p , am ) ; 

cx = 31,0 ; cy = 40.0 ; 
movtcurabs ( acx , acy ) ; 

text ( “SEARCH RADAR PARAMETERS") ; 
settextcl r ( apurple , ablack ) ; 

cx= 5.0 ; cy = 35.0; 
movtcurabs ( acx , acy ) ; 

text("For the explanations of these parameters please 

cy = 32.0 ; 

movtcurabs ( acx , acy ) ; 

text( "refer to the manual . “ ) ; 

h = 1 ; 

settext ( ah , aw , ap , am ) r 
cx = 5.0 ; cy = 23.0 ; 

movtcurabs ( acx , acy ) ; 



text ( "a . PJ = 


k . MF = " ) 


cy = 21 . 0 ; 




movtcurabs ( acx , acy ) ; 




text ( "b. RT 


l.PL = ") 


cy = 19.0 ; 




movtcurabs ( acx , acy ) ; 




text( "c . BAZ = 


m. NS = ") 


cy = 17.0 ; 




movtcurabs ( acx , acy ) ; 




tex t ( "d . GJJ = 


n.BL = " ) 



cy = 15.0 ; 

movtcurabs ( acx , acy ) ; 
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text ( "e . T = 

cy = 13 . 0 ; 

movtcurabs ( Sex , Scy ) ; 

text ( "f . BEL 
cy = 11.0 ; 

movtcurabs ( Sex . Scy ) ; 

text ( "g . JR = 

cy = 9.0 ; 

movtcurabs ( Sex , Scy ) ; 

text ( "h . AZR = 
cy = 7.0 ; 

movtcurabs ( Sex , Scy ) ; 

text ( " i . POWER = 
cy = 5.0 ; 

movtcurabs ( Sex , &cy ) ; 

text ( " j . FAN = 
cx = 45.0 ; cy = 3.0 ; 
movtcurabs ( Sex , Scy ) ; 

text ( "x . EXIT " ) ; 
pointerl = f open ( "search . df 1 
f scant ( poin ter 1 , “%f %f %f %f 
%f %f %f %f 
%f " , 



o . EA = " ) ; 

P.AA = ") ; 

q. ST = ") ; 

r . PRF = " ) ; 

s . BW = " ) ; 

t. FF = ") ; 

, "r " ) ; 

%f %f %f %f %f %f %f %f %f 



%f %f 



&PRF . &AA , ST , &RT , SEA , SFF , SPL . SBL , SBW , SST . SNF , SPOWEP. . &AZE , &NS , SBAZ , SBE 
SFAN , SP J , SGJJ , & JR ) ; 

sprintf ( s , "%f " , PJ) ; 
cx = 21.0 ; cy = 23.0 ; 
movtcurabs ( Sex , Scy ) ; 

text(s) ; 

sprintf (s , "%f" , NF ) ; 

cx = 65.0 ; 
movtcurabs ( Sex , Scy ) ; 

text ( s ) ; . 



sprintf ( s , 


"%f" , 


RT ) ; 


cx = 21.0 ; 


cy =21 


.0 ; 


movtcurabs ( Sex , Scy ) 


i 


text(s) ; 

sprintf ( s , 
cx = 65.0 ; 


"%f" , 


PL ) ; 


movtcurabs ( Sex , Scv ) 


» 


text(s) ; 
spr intf ( s , 


"%f" , 


BAZ ) 


cx = 21.0 ; 


cy = 1 9 


. 0 ; 


movtcurabs ( Sex , Scy ) 


i 


text(s) ; 

sprintf ( s , 
cx = 65.0 ; 


"%f" , 


NS ) ; 


movtcurabs (Sex, Scy) 


f 


text ( s ) ; 

spr intf ( s , 




GJJ ) ; 


cx = 21.0 ; 


cy = 17 


.0 ; 
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movtcurabs ( Sex , Scy ) ; 

text ( s ) ; 



sprintf ( s , 
cx = 65,0 ; 




BL ) ; 


movtcurabs ( & cx, &cy ) 


i 


text (s ) ; 

sprintf ( s , 


'■%f , 


T ) ; 


cx = 21.0 ; 


cy = 1 5 


.0 ; 


movtcurabs! Sex . Scy ) 


i 


text ( s ) ; 

spr intf ( s , 
cx = 65.0 ; 


'0 t 


EA ) ; 


mov tcu rabs ( 


Sex , Scy ) 




tex t ( s ) ; 

spr intf ( s , 


'•%f" , 


BEL ) ; 


cx = 21.0 ; 


cy = 1 3 


.0 ; 


movtcurabs ( Sex , Scy ) 


i 


text ( s ) ; 

sprintf ( s , 
cx = 65.0 ; 




AA ) ; 


movtcurabs ( Sex , Sc:/ ) 


i 


text(s) ; 
sprintf ( s , 


'•%f' , 


JE ) ; 


cx = 21.0 ; 


cy = 1 1 


.0 ; 


mov tcu rabs ( Sex , Scy ) 


i 


text(s) ; 
sprin tf ( s , 
cx = 65.0 ; 


'* 9-f " 

O -L i 


ST ) ; 


movtcurabs ( Sex , Scy ) 




text(s) ; 
spr intf ( s , 


"*£“ , 


AZR ) ; 


cx = 21.0 ; 


cy = 9. 


0 ; 


movtcurabs ( Sex , Scy ) 


i 


text ( s ) ; 

spr intf ( s , 
cx = 65.0 ; 


'•%f" , 


FRF ) ; 


movtcurabs ( Sex , Scy ) 


t 


text(s) ; 
spr intf ( s , 


•'%f" , 


POWER ) 


cx = 21.0 : 


cy = 7 . 


0 ; 


movtcurabs (Sex, Scy) 




text(s) ; 
spr in tf ( s , 


"%f M , 


BW ) ; 



cx = 65.0 ; 
movtcurabs ( Sex , Scy ) ; 

tex t ( s ) ; 

sprintf(s ,"%f" , FAN ) ; 

cx = 21.0 ; cy =5.0 ; 
movtcurabs! Sex , Scy ) ; 

text ( s ) ; 

sprintf (s , "%f" , FF ) ; 
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cx = 65.0 ; 
movtcurabs ( &cx , &cy ) ; 

text(s) ; 
return ; 
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changel ( ) 

{ 

int c ; 

FILE *pointer2 ; 

static char filename[13] ; 

h = 2 ; w = 1 ; 
p ’= 0 ; m = 0 ; 
settext ( &h , &w , &p , &m ) ; 

cx = 5,0 ; cy = 29.0 ; 
movtcurabs ( Sex , &cy ) ; 

text( "Would you like to change any of these values ? 
") ; ' 
cx = 93.0 ; cy = 29.0 ; 
movtcurabs ( &cx , &cy ) ; 

i = getch( ) ; 

if ( i == ‘n‘ ) 

text ( "No " ) ; 

if ( i == 'y' ) 

i 

text ( "Yes " ) ; 

while ( i == ’ y' ) 

i 

h = 1 ; 

settext ( &h , &w , &p , &m ) ; 

c = getch ( ) ; 

swi tch ( c ) 

{ 

case ' a ’ : 

cx = 21.0 ; cy = 23.0 ; 
movtcurabs ( &cx , &cy ) ; 

sprintf(s ,"%f" , PJ ) ; 

text(s) ; 

movtcurabs ( &cx , &cy ) ; 

scanf ( "%f " , &PJ ) ; 

spr in tf ( s , " %f " , PJ ) ; 

text(s) ; 
break ; 
case 'b' : 

cx = 21 . 0 ; cy = 21 . 0 ; 

movtcurabs ( Sex , &cy ) ; 

sprintf(s ,"%f" , RT ) ; 

tex t ( s ) ; 

mov tcurabs ( Sex , &cy ) ; 

scanf ( "%f " , XRT ) ; 

sprintf (s , "%f " , RT ) ; 

text(s) ; 
break ; 
case 'c' : 

cx = 21.0 ; cy = 13.0 ; 
movtcurabs ( &cx , 3-cy ) ; 



<Y/N> 
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sprintf(s ,"%f" , BAZ ) 

text(s) ; 

movtcurabs(8cx, Jcy ) ; 

scanf("%f" , &BAZ ) ; 
sprintf (s , "%f " , BAZ) ; 
text ( s ) ; 

break ; 
case ' d ' : 

cx = 21.0 ; cy = 17.0 ; 

movtcurabs ( gcx , Scy ) ; 

sprintf (s , "%f " , GJJ ) 

text(s) ; 

movtcurabs ( Sex , Scy ) ; 

scanf("%f'' , &GJJ ) ; 
sprintf (s, "%f ", GJJ) ; 
tex t ( s ) ; 

break ; 
case 'e* : 

cx = 21.0 ; cy = 15.0 ; 
movtcurabs ( Sex , Scy ) ; 

sprintf ( s , "%f " , T ) ; 

text ( s ) ; 

movtcurabs ( sex , &cy ) ; 

scanf ( "%f " , &T ) ; 

sprintf ( s , "%f " , T ) ; 

tex t ( s ) ; 

break ; 
case ' f ' : 

cx = 21.0 ; cy = 13.0 ; 
mov tcurabs ( Sex , &cy ) ; 

sprintf (s ,'‘%f" , BEL ) 
text ( s ) ; 

movtcurabs ( Sex , Scy ) ; 

scanf ("%f" , &BEL ) ; 

sprintf (s , "%f " , BEL ) ; 

text(s) ; 
break ; 
case ’ g ’ : 

cx = 21.0 ; cy = 11.0 
movtcurabsUcx.&cy ) ; 

sprintf (s ,"%f" . JR ) 

text(s) ; 

movtcurabs ( Sex , &cy ) ; 

scanf ( " %f " , S JR ) ; 

sprintf (s , "%f" , JR ) ; 

text ( s ) ; 

break ; 
case * h ’ : 

cx = 21.0 ; cy = 3.0 ; 
mov tcu rabs ( Sex , Scy ) ; 

sprintf (s ,"%f" , AZR 



) 



text(s) ; 

mov tcurabs ( &cx , Scy ) ; 
scanf ("%£" , &AZR ) ; 

sprintf ( s , "%f " , AZR ) ; 

text(s) ; 
break ; 
case ' i ’ : 

cx = 21.0 ; cy = 7.0 ; 
movtcurabs ( Sex , &cy ) ; 

sprintf (s ,"%f" , POWER ) ; 

text(s) ; 

movtcurabs ( Sex , &cy ) ; 

scanf (••%£■• , & POWER ) ; 

sprintf (s, "%f", POWER ) ; 

text(s) ; 
break ; 
case ' j ’ : 

cx = 21.0 ; cy = 5.0 ; 
movtcurabs ( &cx , &cy ) ; 

sprintf (s ,"%f" , FAN ) ; 

text(s) ; 

movtcurabs ( &cx , &cy ) ; 

scanf ("%f" , &FAN ) ; 

sprintf (s , "%f " , FAN ) ; 

text(s) ; 
break ; 
case ' k ' : 

cx = 65,0 ; cy =23,0 ; 
movtcurabs ( &cx , &cy ) ; 

sprintf (s ,''%f" , NF ) ; 

text(s) ; 

movtcurabs ( &cx , &cy ) ; 

scanf ("%f" , &NF ) ; 

sprintf (s, "%f ", NF ) ; 

text(s) ; 
break ; 
case ' 1 ' : 

cx = 65.0 ; cy =21,0 ; 
movtcurabs ( &cx , &cy ) ; 

sprintf (s ,"%f" , PL ) ; 

text ( s ) ; 

movtcurabsi Sex , ?<cy ) ; 

scanf ( "%f " , SPL ) ; 

sprintf (s , "%f M , PL ) ; 

text(s) ; 
break ; 
case 'm' : 

cx = 65.0 ; cy = 19,0 ; 
mov tcur abs ( Sex , &cy ) ; 

sprintf (s ,"%f" , NS ) ; 

text ( s ) ; 
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movtcurabs ( gcx , gey ) ; 

scanf ( ”%f " , gNS ) ; 

sprintf ( s , " %f " , NS ) ; 

text(s) ; 
break ; 
case ' n ' : 

cx = 65.0 ; cy = 17.0 
movtcurabs ( gcx , gey ) ; 

sprintf (s ,"%f" , BL ) 
text ( s ) ; 

movtcurabs ( gcx , gey ) ; 

scanf ("%f" , SEL ) ; 

sprintf (s . "%f " . BL ) ; 

text ( s ) ; 

break ; 
case 'o' : 

cx = 65,0 ; cy = 15.0 
mov tcurabs ( gcx , &cy ) ; 

sprintf (s , "%f" , EA ) 
tex t ( s ) ; 

movtcurabs ( gcx , gey ) ; 

scanf ( "%f " , gEA ) ; 

sprintf (s, "%f ", EA ) ; 

text(s) ; 
break ; 
case ' p ' : 

cx = 65.0 ; cy = 13.0 
movtcurabs ( gcx , gey ) ; 

sprintf (s ,"%f" , A A ) 

text(s) ; 

mov tcurabs ( gcx , gey ) ; 

scanf ("sfef" , gAA ) ; 

sprintf ( s , " %f " , AA ) ; 

tex t ( s ) ; 

break ; 
case ' q ' : 

cx = 65.0 ; cy = 11.0 
movtcurabs ( gcx , gey ) ; 

sprintf (s ,"%f" , ST ) 
text ( s ) ; 

movtcurabs ( gcx , gey ) ; 

scanf ( "%f " , gST ) ; 

spr in tf ( s , " %f " , ST ) ; 

text(s) ; 
break ; 
case ’ r ’ : 

cx = 65,0 ; cy =9.0 ; 
movtcurabs ( gcx , gey ) ; 

sprintf (s ,‘'%f" , PRF) 

text ( s ) ; 

movtcurabs ( gcx , gey ) ; 



scanf ( "%f " , &FRF ) ; 

sprintf ( s , "%f " , PRF) ; 
text ( s ) ; 

break ; 
case 's’ : 

cx = 65.0 ; cy = 7.0 ; 
movtcurabsUcx , acy ) ; 

sprintf ( s ,"%f" , BW ) ; 

text ( s ) ; 

movtcurabs ( Sex , Scy ) ; 

scanf ( "%f " , SBW ) ; 

sprintf (s , "%f " , BW ) ; 

text ( s ) ; 

break ; 
case ’ t ' : 

cx = 65.0 ; cy = 5.0 ; 
movtcurabs ( Sex , Scy ) ; 

sprintf (s ,"%f" , FF ) ; 

text(s) ; 

movtcurabs ( Sex , Scy ) ; 

scanf ("%f" , &FF ) : 

sprintf ( s , "%f " , FF ) ; 

text ( s ) ; 

break ; 
case ' x ' : 

h = 2 ; w = 1 ; 

settext ( Sh , &w , &p , &m ) ; 

cx = 93 . 0 ; cy = 29 . 0 ; 
movtcurabs ( Sex , Scy ) ; 

text ( "Yes " ) ; 

cx = 5.0 ; cy = 29.0 ; 
movtcurabs ( Sex , Scy ) ; 

text( "Would you like to change any of these values ? <Y/N> “ ) 



movtcurabs ( Sex , Scy ) ; 

text ("DO YOU WISH TO SAVE YOUR DATA YOU JUST ENTERED ?") ; 

i = getch (); 
if ( i == ' y ' ) 

movtcurabs ( Sex , key ) ; 

tex t ( "DO YOU WISH TO SAVE YOUR DATA YOU JUST ENTERED ?") : 



mov tcu rabs ( Sex , Scy ) ; 

tex t ( "WHAT IS YOUR DATA FILENAME ? <f i 1 ename . ex t> ") ; 

scanf("%s", filename) ; 
text ( f i 1 ename ) ; 



"w " ) 
%f %f 



%f %f %f %f %f %f %f %f 



PRF 



pointer2 = f open ( f i 1 ename 
f printf ( pointer2 , "%f %f %f 
%f %f %f %f %f %f", 

AA , T . RT , EA , FF . PL , BL , BW , ST , IMF , POWER , AZR , NS , BAZ , BEL , FAN , PJ , GJJ 



JR) 



i 



l = 



n 
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} 

} 

} 

h = 1 ; 

settext ( &h , &w , &p , jm ) ; 

cx = 60.0 ; 

cy = 1.0 ; 

settextcl r ( &yel 1 ow , &black ) ; 

movtcurabs ( &cx , &cy ) ; 

text(''Hit any Key to continue") ; 

i = getch ( ) ; 

setcol or ( Sblack ) ; 

i = getch ( 12 ) ; 

cl r ( ) ; 

return ; 

1 
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input2 ( ) 

{ 

FILE *pointer3 ; 



setviewport(svxl , S vy 1, Svx2, Svy2,S bo rder.s back) ; 

pixels = 13 ; 

set 1 nwidth ( Spixe 1 s ) ; 

setco 1 or ( Sred ) ; 

box ( Sxl , Syl , Sx2 , Sy2 ) ; 

cx = 50.0 ; 

cy = 49 . 0 ; 

movabs ( Sex , Scy ) ; 

f 1 ood ( & 1 ibl ue ) ; 

setxor ( &swi tsh ) ; 

h = 3 ; w = 2 ; 

p = 0 ; m = 0 ; 

settext ( Sh , Sw , Sp , Sm ) ; 

settextc 1 r ( Sgreen , Sbl acK ) ; 

cx = 30.5 ; cy = 45.0 ; 

movtcurabs ( Sex , Scy ) ; 

text ("INPUT MODULE”) ; 

settextcl r ( Sbl ue , Sblack ) ; 

h = 2 ; M = 1 ; 

set tex t ( Sh , Sw , S p , S m ) ; 

cx = 30.0 ; cy = 40,0 ; 

movtcurabs ( Sex , Scy ) ; 

text ( "TRACKING RADAR PARAMETERS") ; 

settextcl r ( Spurple , Sblack ) ; 

cx = 5.0 ; cy = 35,0; 

movtcurabs ( Sex , Scy ) ; 



text("For the explanations of 
cy = 32. 0 ; 
movtcurabs ( Sex , Scy ) ; 
text( "refer to the manual. ") 

h = 1 ; 

settext ( Sh , S’w , Sp , Sm ) ; 

cx= 5.0 ; cy =23.0; 
mov tcurabs ( Sex , Scy ) ; 


these parameters 
» 


text ("a. EAT 

cy = 2 1 . 0 ; 

mov tcurabs ( Sex , Scy ) ; 


g.CST = 


" ) ; 


text ( "b . TAZ 
cy = 1 9 . 0 ; 

movtcurabs ( Sex , Scy ) ; 


h . NFT = 


" ) ; 


text ( ”c . GLT = 

cy = 17.0 ; 

movtcurabs ( Sex , Scy ) ; 


i.BLT = 


" ) ; 


text ( "d . TBEL = 
cy = 15.0 ; 


j . AAT = 


" ) ; 
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X. BWT 



movtcura bs ( Sex , Scy ) ; 

tex t ( "e . STT = 

cy = 1 3 . 0 ; 

movtcura bs ( Sex , Scy ) ; 

tex t ( "f . POWT = 
cy = 11.0 ; 
movtcurabs ( Sex , Scy ) ; 

text ( "x. EXIT" ) ; 



1 . FFT = " ) 



pointer3 = f open ( "track . df 1 " , " r ” ) ; 

f scant ( pointer3 , ”%f %f %f %f %f %f %f 

%f " , 4 A AT . SEAT . 3 BWT , 3CST , 

&FOWT , &FFT . SGLT , &BLT . S3TT , SNFT , 3TAZ , ST BEL ) 

sprintf ( s , "%f " , EAT) ; 

cx = 21.0 ; cy = 23.0 ; 

movtcurabs ( &cx , Scy ) ; 

text(s) ; 

sprintf (s , , CST ) ; 

cx = 65.0 ; 
movtcurabs ( Sex , Scy ) ; 

text (s) ; 

sprintf ( s , "%f" , TAZ ) ; 

cx = 21.0 ; cy = 21.0 ; 
movtcurabs ( Sex , Scy ) ; 

text(s) 

sprintf (s , "%f" , NFT ) ; 

cx = 65.0 ; 
movtcurabs ( Sex , Scy ) ; 

text ( s ) ; 

sprintf (s , "%f " , GLT ) ; 

cx = 21.0 ; cy = 19.0 ; 
movtcurabs ( Sex , Scy ) ; 

text ( s ) ; 

sprintf (s , "%f" , BLT ) ; 

cx = 65 . 0 ; 
movtcurabs ( Sex , Scy ) ; 

text(s) ; 

sprintf (s , "%f" ,TBEL) ; 

cx =21.0 ; cy = 17.0 ; 
movtcurabs ( Sex , Scy ) ; 

text ( s ) ; 

sprintf (s , "%f" , AAT ) ; 

cx = 65.0 ; 
movtcurabs ( Sex , Scy ) : 

text ( s ) ; 

sprintf ( s , "%f" , STT) ; 

cx =21.0 ; cy = 15.0 ; 
movtcurabs ( Sex , Scy ) ; 

text ( s ) ; 

sprintf (s , "%f” , BWT); 

cx = 65.0 ; 
movtcurabs ( Sex , Scy ) ; 



%f 



%f 



2-f 

O 
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text ( s ) ; 

sprintf(s , "%f" , POWT ) ; 

cx = 2 1 . 0 ; cy = 1 3 . 0 ; 

movtcurabs ( &cx , &cy ) ; 

text(s) ; 

sprintf(s , ”%f" , FFT ) ; 

cx = 65.0 ; 
movtcurabsf Sex , &cy ) ; 

text(s) ; 
return ; 
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change2 ( ) 

{ 

int c ; 

static char filename[13] ; 

FILE *pointer4 ; 

h = 2 ; 
w = 1 ; 

P = 0 ; 

m = 0 ; 

settext ( &h , Sw , &p , Sm ) ; 

cx = 5.0 ; cy = 29.0 ; 
movtcurabs ( Sex , &cy ) ; 

text( "Would you like to change any of these values ? <Y/N> 
» 

cx = 93.0 ; cy = 29.0 ; 
movtcurabs ( Sex , Scy ) ; 

i = getchl ) ; 

if ( i == 'n' ) 

text ( "Mo " ) ; 

if ( i = = ' y ' ) 

s 

< 

text ( " Yes " ) ; 

while ( i == 'y' ) 

{ 

h = 1 ; 

settext ( Sh , Sw , Sp , Sm ) ; 

c = getch ( ) ; 

switch ( c ) 

case ' a ' : 

cx = 21,0 ; cy = 23.0 ; 
movtcurabs ( Sex , Scy ) ; 

sprintf(s ,"%f" , EAT) ; 
text ( s ) ; 

movtcurabs ( Sex , &cy ) ; 

scanf("%f" , SEAT) ; 
sprintf ( s , "%f " , EAT) ; 
text ( s ) ; 

break ; 
case ‘ b 1 : 

cx = 21.0 ; cy = 21.0 ; 

movtcurabs ( Sex , Scy ) ; 

sprintf (s ,"%f" , TAZ) ; 
text ( s ) ; 

mov tcurabs ( Sex , Scy ) ; 

scanf("%f" , STAZ) ; 
sprintf ( s , "%f " , TAZ) ; 
text ( s ) ; 

break ; 
case 'o' : 
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cx = 21,0 ; cy = 19.0 ; 
movtcurabs ( &cx , &cy ) ; 

sprintf(s ,"%f" , GLT ) 
text(s) ; 

mov tcurabs ( &cx , &cy ) ; 

scanf("%f“ , &GLT ) ; 

sprintf ( s , "%f " , GLT) ; 
tex t ( s ) ; 

break ; 
case 'd' 

cx = 21.0 ; cy = 17.0 ; 

mov tcu rabs ( .Sex , &cy ) ; 

sprintf (s , "%f" , TBEL ) 
text(s) ; 

movtcurabs ( & c x , &cy ) ; 

scanf ( "%f " , &TBEL ) ; 

sprintf (s , "%f ", TBEL) ; 
text ( s ) ; 

break ; 
case 'e’ : 

cx =21,0 ; cy = 15.0 ; 

movtcurabs ( & c x , &cy ) ; 

sprintf (s , "%f" ,STT) ; 
text ( s ) : 

mov tcu rabs ( &cx , &cy ) ; 

scanf ("%f" , &STT ) ; 

sprintf (s , "%f ", STT ) ; 

text(s) ; 
break ; 
case ' f ' : 

cx = 21.0 ; cy = 13.0 ; 
movtcurabs ( &cx , &cy ) ; 

sprintf (s , "%f" , FOWT ) 
text(s) ; 

movtcurabs ( &cx , &cy ) ; 

scanf ("%f" , &POWT ) ; 

sprintf ( s , "%f " , POWT ) ; 

text(s) ; 
break ; 
case ’ g ' : 

cx = 65.0 ; cy = 23.0 ; 
mov tcurabs ( &cx , 3cy ) ; 

sprintf (s ,“%f" , GST ) 

text ( s ) ; 

mov tcu rabs ( &cx , &cy ) ; 

scanf ("%f" , &CST ) ; 

sprintf (s , "%f " , CST ) ; 

.text(s) ; 
break ; 
case ' h ' : 

cx = 65.0 ; cy = 



21 . 0 



movtcurabs ( Sex , Scy ) ; 

sprintf(s , "%f" , NFT ) ; 

text(s) ; 

movtcurabs ( Sex , Scy ) ; 

scanf ( " % f " , SNFT ) ; 

sprintf ( s , "%f " , NFT ) ; 

text ( s ) ; 

break ; 
case ' i ' : 

cx = 65,0 ; cy = 19.0 ; 
movtcurabs ( Sex , Icy ) ; 

sprintf (s ,"%f" , ELT ) ; 

text ( s ) ; 

movtcurabs ( Sex , Scy ) ; 

scanf ("%f" , SBLT ) ; 
sprintf ( s , "%f " , BLT ) ; 

text(s) ; 
break ; 
case ' j ‘ : 

cx = 65.0 ; cy = 17.0 ; 
movtcurabs ( Sex , Scy ) ; 

sprintf (s ,"%f" , AAT ) ; 

text(s) ; 

movtcurabs ( Sex , Scy ) ; 

scanf ("%f" , SAAT ) ; 
sprintf (s , "%f " , AAT ) ; 

text ( s ) ; 

break ; 
case ' k ' : 

cx = 65.0 ; cy = 15.0 ; 
movtcurabs ( Sex , Scy ) ; 

sprintf (s ,"%f" , BWT) ; 
text(s) ; 

movtcurabs ( Sex , Scy ) ; 

scanf ("%f'‘ , SBWT ) •, 

sprintf ( s , "%f " , BWT) ; 

text(s) 

break ; 

case * 1 * : 

cx = 65.0 ; cy = 13.0 ; 
movtcurabs ( Sex , Scy ) ; 

sprintf (s , "%f" , FFT ) ; 

tex t ( s ) ; 

movtcurabs ( Sex , Scy ) ; 

scanf ("%f" , SFFT) ; 

sprintf ( s , "%f " , FFT ) ; 

text ( s ) ; 

break ; 

case ' x ’ : 

h = 2 ; w = 1 ; 

set text ( Sh , Sw , Sp‘, Sm ) ; 
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cx = 10.0 ; cy = 8.0 ; 

movtcurabs ( Jcx , &cy ) ; 

text ("DO YOU WISH TO SAVE YOUR DATA YOU JUST ENTERED ?") ; 

i = getch ( ) ; 

if (i == 'y ■ ) 

{ 

mov tcurabs ( &cx , &cy ) ; 

text ("DO YOU WISH TO SAVE YOUR DATA YOU JUST ENTERED ?") 
movtcurabs ( 8-cx , &cy ) ; 

text ("WHAT IS YOUR DATA FILENAME ? <f i 1 ename . ext> ") : 

scanf("%s", filename) ; 
text ( f i 1 ename ) ; 

pointer4 = f open ( f i 1 ename , "u") ; 

f printf ( pointer4 , ”%f %f %f %f %f %f %f %f %f %f %f 
%f " , AAT , EAT . BWT , 

CST , POWT , FFT , GLT . BLT , STT , NFT , TAZ . TBEL ) ; 

} 

i = ' n ' ; 

} 

} 

1 

r 

h = 1 ; 

settext ( &h , &w , &p , &m ) ; 

cx = 60.0 ; 

cy = 1.0 ; 

settextcl r ( &yel 1 ow , ablack ) ; 

movtcurabs ( &cx , &cy ) ; 

text("Hit any key to continue") ; 

i = getch ( ) ; 

setcol or ( &bl ack ) ; 

cl r ( ) ; 

return ; 
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^include "stdio.h 
^include "math.h" 



float xl = 0.0 , yl = 0.0 
float x2 = 1 00 . 0 , y2 = 100.0 ; 

float v 1 1 = 0.635 . vl2 = 0.0 , vl3 
float v2 1 = 0.0 . v22 = 0.0 , v23 = 
float max_range , lrange , mrange , 
float dlrange, dmrange , dsrange ; 



float cx 



cy 



int 


back = 


-1 


, b 


int 


h , W , 


p 


' , m 


int 


swi tsh 


= 


1 ; 


int 


black 


— 


0 ; 


int 


dbl ue 


= 


1 ; 


int 


green 


= 


2 ; 


int 


1 bl ue 


= 


3 ; 


int 


red 


= 


4 ; 


int 


purple 


= 


5 ; 


int 


brown 


= 


6 ; 


int 


gray 




7 ; 


int 


dgray 


= 


8 ; 


int 


blue 


= 


9 : 


int 


1 green 


= 


10 ; 


int 


v 1 bl ue 


= 


11 ; 


int 


pink 


= 


12 ; 


int 


lpurple 




13 ; 


int 


yellow 


- 


14 ; 


int 


white 


= 


15 ; 


float RCPA 


i 


RB ; 



= 1.0 , 
0.635 . 
srange 



vl 4 
v24 



1.0 ; 

1.0 ; 



float xx , yy 
float XD , YD 
float xt , y t 
float zt 
float U , at 
float R 



as ; /* initial position of the target */ 
/* constant altitude of the target */ 

/* velocity and course of target */ 

/* distance between ship and target*/ 



main ( ) 

{ 

int mode 



setdev ( "HALOIBME. DEV " ) ; 

mode = 4 ; 

ini tgraphics ( Xmode ) ; 

setwor 1 d ( &xl , &y 1 , &x2 , &y2 ) ; 

CPA ( ) ; 

Input ( ) ; 

Opening ( ) ; 

inf o ( ) ; 
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draw ( ) ; 

detect ( ) ; 

c 1 osegraphics ( ) ; 






} 






Inpu t ( ) 

{ 



printf ( "\n PLEASE ENTER THE MAXIMUM RANGE DETECTION ” ) ; 

scanf("%f" , &max_range ) ; 

printf ( "\n PLEASE ENTER THE LONG , MEDIUM AND SHORT RANGE 
DISTANCES \n" ) ; 

scanf ( "%f %f %f " , slrange , Smrange , asrange ) ; 

dlrange = 45.0 * lrange / max_range ; 

dmrange = 45.0 * mrange / max_range ; 

dsrange = 4.5.0 * srange / max_range ; 

se tco 1 or ( &bl ack ) ; 

clr() ; 

return ; 



i 
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Opening ( ) 

* 



float bx , by ; 

int j ; 

int pixels ; 

float radius ; 

float aspect = 0.365 ; 

pixels=13; 
set 1 nwidth ( apixel s ) ; 

setcol or ( ared) ; 
box ( axl , ayl , &x2 , &y2 ) ; 

setviewport(avll,avl2,avl3,avl4,aborder, aback ) 
setxor ( aswi tsh ) ; 

pixels = 1 ; 

set 1 nwidth ( apixels ) ; 

cx = 0 . 0 ; cy = 0 . 0 ; 
movabs ( acx , acy ) ; 

cy = 100.0 ; 

1 nabs ( a cx , acy ) ; 

cx = 50.0 

cy = 98 , 0 ; 

movabs ( acx , acy ) ; 

flood (alblue) ; 

setvieuiport( av21 ,av22 ,av23 ,av24 , aborder , aback ) 
setxor ( aswi tsh ) ; 

cx = 49,0 ; cy = 49,0 ; 
bx = 51.0 ; by = 51.0 ; 
box ( acx . acy , abx , aby ) ; 

setcol or ( ayel 1 ow ) ; 

cx ='1,0 ; cy = 50.0 ; 
movabs ( acx , acy ) ; 

cx = 99 . 0 ; 

1 nabs ( acx , acy ) ; 

cx = 50.0 ; cy = 98.5 ; 

movabs ( acx , acy ) ; 

cy = 1 . 5 ; 

lnabs ( acx , acy ) ; 

cy = 0.0; 

f or ( j =0 ; j <19 ; j ++ ) 

cx = 49 . 5 ; 

cy = cy +5.0 ; 
movabs ( acx , acy ) ; 

cx =50.5 ; 

1 nabs ( acx , acy ) ; 

} 

cx = 0.0 ; 

fort j=0;j<19;j++) 

{ 
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cx = cx + 5 . 0 ; 
cy = 49.5 ; 

movabs ( Sex , &cy ) ; 

cy = 50 . 5 ; 

1 nabs ( Sex , &cy ) ; 

} 

setcolor(Sgreen) ; 

cx = 50.0 ; cy = 50.0 ; 

movabs ( Sex , &cy ) ; 

setasp ( ^aspect ) ; 

radius =45.0 ; 

ci r ( Sradius ) ; 

setcol or ( Sdblue ) ; 

movabs ( Sex , &cy ) ; 

ci r ( &dl range ) ; 

movabs ( &cx , &cy ) ; 

cir(Xdmrange) ; 

movabs ( Sex , &cy ) ; 

ci r ( Sdsrange ) ; 

h = 1 ; m = 1 ; 

p = 0 ; m = 0 ; 

set tex t ( &h , &w , &p , &m ) ; 

se ttex tc 1 r ( Sgreen , &bl acK ) ; 

cx = 83.0 ; cy = 83.0 ; 

movtcurabs ( &cx , &cy ) ; 

text ( "DR " ) ; 

cx = 50.0 + dlrange ; 

cy = 51.0 ; 

mov tcurabs ( &cx , &cy ) ; 

text ( "LR " ) ; 

cx = 50,0 + dmrange ; 

movtcurabs ( &cx , Scy ) ; 

text ( "MR" ) ; 

cx = 50.0 + dsrange ; 

movtcurabs ( &cx , &cy ) ; 

text ( "SR" ) ; 

setviewport(&vll,&vl2,&vl3,&vl4, & border , &back ) 
setxor ( Sswi tsh ) ; 

setcolor(Sbrown) ; 
cx = 0.0 ; cy = 70.0 ; 
movabs ( Sex , Scy ) ; 

cx = 100.0 ; 
lnabs( Sex , &cy ) ; 

cy = 40 ; 
movabs ( Sex , &cy ) ; 

cx = 0.0 ; 

1 nabs( &cx , &cy ) ; 

h = 1 ; w •= 1 ; 
p = 0 ; m = 0 ; 
settext ( &h , &w , &p . Sm ) ; 
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settextcl r ( Slgreen , Sblack ) ; 

cx = 40.0 ; cy = 95.0 ; 

movtcur abs ( Sex , scy ) ; 

text ( "RADAR " ) ; 

cx = 35.0 ; cy = 67.0 ; 

movtcurabs ( Sex , Scy ) ; 

text ( "MISSILES" ) ; 

cx = 38.0 ; cy = 37.0 ; 

movtcurabs ( Sex , Scy ) ; 

text( "TARGET” ) ; 

settextcl r ( Sdgray , Sblack ) ; 

cx = 13.0 ; cy = 90,0 : 

mov tcurabs ( Sex , Scy ) ; 

text( "Search Radar active"); 

cx = 5 . 0 ; cy = 84 . 0 ; 

settextcl r ( Spink , Sbl ac!c ) ; 

movtcurabs ( Sex , Scy ) ; 

text ( "Frequency : MHZ") 

cy = 80.0 ; 

movtcurabs ( Sex , Scy ) ; 

text( "Max detec : " ) ; 

cx = 28.0 : cy = 72.0 ; 

movtcurabs ( Sex , Scy ) ; 

settextcl r ( Sgray , Sbl ack ) ; 

text ( " MO CONTACT ") ; 

cx = 5.0 ; cy = 63.0 ; 

movtcurabs ( Sex , Scy ) ; 

settextcl r ( Spink , Sblack ) ; 

text ( " RANGE QTY STATUS") 

cx = 5.0 ; cy = 58.0 ; 

movtcurabs ( Sex , Scy ) ; 

text ( "Long : “ ) ; 

cy = 54.0 ; 

movtcurabs ( &cx , scy ) ; 

text ( "Medium ) ; 

cy = 50,0 ; 

movtcurabs ( Sex , Scy ) ; 

text ( "Short : " ) ; 

cy = 32 . 0 ; 

movtcurabs ( Sex , Scy ) ; 

text ( " ID : " ) ; 

cy = 29.0 ; 

movtcurabs ( Sex , Scy ) ; 

text ( " X : " ) ; 

cy = 26.0 ; 

movtcurabs ( Sex , Scy ) ; 

text ( " Y :") ; 

cy = 23,0 ; 

movtcurabs ( Sex , Scy ) ; 

tex t ( " Z : ” ) ; 

cy = 20. 0 ; 
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movtcurabs (Sex, Scy) 
text ("RANGE : " ) ; 

cy = 17.0 ; 
movtcurabs ( Sex , &cy ) 
text ( “ CPA :") ; 

cy = 14.0 ; 
movtcurabs ( Sex , &cy ) 
text ( ” X : ” ) ; 
cy = 11.0 ; 
mov tcurabs ( sex . &cy ) 
text ( " Y ; 

cy = 8.0; 

movtcurabs( Sex , Scy ) 
text ( "Speed ; " ) ; 

cy = 5 . 0 ; 
movtcurabs ( Sex , &cy ) 
text ( "Course : " ) ; 

return ; 



inf o ( ) 

float freq = 15.5679 ; 

int sqty = 16 . mqty = 8 , lqty = 8 ; 

int i ; 

static char s[14] ; 

setvie wport( &v 11, a vl2, & v 13, Svl4,S border, aback) 
setxor ( 8swi tsh ) ; 

settextcl r ( Sgray , Sblack ) ; 

cx =42.0 ; cy = 84.0 ; 
movtcurabs ( &cx , &cy ) ; 

sprintf (s , "%f " , freq) ; 
tex t ( s ) ; 

cy = 80.0 ; 
movtcurabs ( Sex , Xcy ) ; 

sprintf (s , "%f ", max_range ) ; 

text ( s ) ; 

cx = 29.0 ; cy = 58.0 ; 
moutcurabs ( Sex , Scy ) ; 

sprintf (s , "%f ", 1 range ) ; 

text ( s ) ; 

cx = 65.0 ; 
movtcurabs ( &cx , &cy ) ; 

sprintf ( s , "%d" , lqty ) ; 

text(s) ; 
cx = 80.0 ; 
movtcurabs ( &cx , &cy ) ; 

text ( "Ready " ) ; 

cx = 29.0 ; cy = 54.0 ; 
movtcurabs ( Sex , &cy ) ; 

sprintf ( s , "%f *’ .mrange ) ; 

text(s) ; 
cx = 65.0 ; 
movtcurabs ( sex , &cy ) ; 

spr in tf ( s , " %d " , mq ty ) ; 

text(s) ; 

cx = 80 . 0 ; ' ‘ . ■ 

movtcurabs) Sex , Scy ) ; 

tex t ( "Ready " ) ; 

cx = 29.0 ; cy = 50.0 ; 

movtcurabs ( Sex , Scy ) ; 

spr intf ( s , "%f " , srange ) ; 

tex t ( s ) ; 

cx = 65,0 ; 

mov tcu rabs ( Sex , &cy ) ; 

spr intf (s , "%d" , sqty ) ; 

text ( s ) ; 

cx = 80.0 ; 

movtcurabs ( Sex , &cy ) ; 

text ( "Ready " ) ; 
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settextcl r ( &dgray , &bl ac(C ) ; 

cx = 15.0 ; 

cy = 2.0 ; 

movtcurabs ( Sex , &cy ) ; 

text(”Hit any key to continue") ; 

i = getch ( ) ; 

return ; 
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CPA( ) 
i 



float 


xs , 


ys 


/* 


initial position of the ship 


*/ 




float 


vs , 


as 


/* 


velocity and course of the ship 


*/ 




float 


A 




/ * 


angle between horizontal and R 


*/ 




float 


xl . 


y 1 


/* 


define ship’s velocity vector 


*/ 




float 


x2 . 


72 


/* 


define target's velocity vector 


*/ 




float 


V 

■r* t 


y 


/* 


target's relative velocity vect 


*/ 




float 


B 




/* 


angle of x-y vector 


*/ 




float 


C 




/* 


angle between R and x-y vector 


*/ 




float 


CPA 




/'* 


CPA in X-Y plan 






f 1 03 t 


A A 


, 3E 


. CC ; 






rintf ( 


"\n 


PLEASE 


INPUT X, Y AND Z POSITION OF THE TARGET 


scanf ( 


"%f %f %f 1 


* , s 


xt , &y t , Xzt ) ; 







printf ( "\n PLEASE ENTER ITS SPEED AND COURSE 
scanf ( "%f %f " , &U, &at) ; 



") 



printf ("\n PLEASE INPUT X AND Y POSITION OF THE SHIP : " ) 

scanf (" %f %f " , ^xs, Xys) : 

printf ( "\n PLEASE ENTER ITS SPEED AND COURSE : " ) ; 

scanf (" %f %f " , Svs , &as ) ; 



at = at * PI / 130.0 ; 

as = as * PI / 180.0 ; 

R = pow ( ( xs-xt ) * ( xs-xt ) +( ys-y t ) * ( ys-y t ) , 0.5) ; 

A = atan2 ( (y t-ys ) , ( xt-xs ) ) ; 

xl = vs * cos(as) ; 

yl = vs * sin(as) ; 

x2 = U * cos ( at ) ; 

y2=U *sin(at) ' ; 

x = xl - x2 + 0.000001 ; /* to avoid division by zero */ 

y = yl _ y2 ; 

B = atan2 ( y , x ) ; 

C = A - B ; 

if ( C > 1.57079 :i C < -1.57079 ) 

{ 

return ; 

} 

CPA = R * sin ( C ) ; 

ECPA= pow ( ( CPA * CPA ) +( zt * zt), 0.5 ) ; 

AA = PI/2 - A ; 

CC = PI/2 - C ; 

BB = AA - CC ; 

xx = CPA * sin(BB) + xs ; 

yy = CPA * cos(BB) + ys ; 

RB = B + PI ; 

return 
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draw ( ) 

i 

float dxx , dyy ; 

setviewport(&v21,&v22,.Sv23,&v24,& border, a back) ; 
setxor ( &swi tsh ) ; 

XD = 50.0 + 45.0 * xt / max_range ; 

YD = 50.0 + 45.0 * yt / max_range ; 
setco 1 or ( ared ) ; 
movabs ( aXD , aYD ) ; 

dxx = 50.0 + 45.0 * xx / max_range ; 
dyy = 50.0 •+■ 45.0 * yy / max_range ; 

1 nabs ( adxx . adyy ) ; 
cx = 50.0 ; cy = 50.0 ; 
lnabs ( acx , acy ) ; 

return ; 
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detect ( ) 

1 

static char s[14] ; 

inti ; 

float atd ; 

setviewpor t ( &v21 , Sv22 , Sv23 , &v24 , ^border , Sback ) 
setco 1 or < Swhi te ) ; 

while(R > max_range ) 

i 

xt = xt + 10.0 / 3600.0 * U * cos ( RB ) ; 

yt = yt + 10,0 / 3600.0 * U * sin(RB) ; 

R = pow(xt * xt yt * yt + zt * zt , 0.5) 
XD = 50.0 + 45.0 * xt / max_range ; 

YD = 50.0 + 45.0 * yt / max_range ; 

ptabs ( &XD , &YD ) ; 

} 

setviewpor t ( &vl 1 , Svl2 , Svl3 , Svl4 , Sborder , &back ) 
settextcl r( &gray , &black ) ; 

cx = 28.0 ; cy = 72.0 ; 
movtcurabs( Sex , Scy ) ; 

text ( " NO CONTACT") ; 
m = 1 ; 

settext ( sh , &w , sp , &m ) ; 

settextclr(Sl green , &gr ay ) ; 
cx = 1 8 . 0 ; 
movtcurabs ( &cx , Scy ) ; 

text ( "WARNING ! CONTACT") 
m = 0 ; 

settext ( Sh , &w , &p , &m ) ; 

settextcl r ( Sgray , &black ) ; 

cx = 35.0 ; cy = 32.0 ; 
movtcurabst Sex , &cy ) ; 

text ("Hostile") ; 
cy = 17.0 ; 
movtcurabs ( &cx , &cy ) ; 

sprintf (s, "%f " , RCPA ) ; 

text ( s ) ; 

cy = 14.0 ; 
movtcurabs( Sex , Scy ) ; 

sprintf ( s , "%f ", xx ) ; 

text(s) ; 
cy = 11.0 ; 
movtcurabs ( Sex , Scy ) ; 

sprintf ( s , "%f ", yy ) ; 

text(s) ; 
cy = 8.0 ; 

movtcurabs ( Sex , Scy ) ; 

sprintf (s, "%f" ,U) ; 

text ( s ) ; 

cy = 5.0 ; 
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movtcurabs( Sex , Scy ) ; 

atd = 180.0 * at / FI ; 
sprintf ( s , "%f " , atd ) ; 

text(s) ; 

set textc 1 r ( & 1 green , Sbl ack ) ; 

cy = 29.0 ; 

movtcurabsUcx , Scy ) ; 

sprin tf (s , "%f " , xt ) ; 

text ( s ) ; 

cy = 26.0 ; 

movtcurabs(&cx,lcy) ; 

sprintf (s , "%f " ,vt ) ; 

text ( s ) ; 

cy = 23.0 ; 

mov tcu rabs ( Sex , Scy ) ; 

sprintf ( s , "%f ", zt ) ; 

text(s) ; 

cy = 20.0 ; 

movtcu rabs ( &cx , Scy ) ; 

sprintf (s , "%f ", R ) ; 

tex t ( s ) ; 

while(R <= max_rar.ge ) 

{ 

set viewport ( &v21 , &v22 , Sv2 3 , &v24 , Sborder , &back ) 
setcol or ( Swhi te ) ; 

ptabs ( &XD , &YD ) ; 

setviewport ( & vll ,&vl2 ,&vl3 , Svl4 , ^border , &back ) 
setxor ( Sswi tsh ) ; 

m = '1 ; 

settext ( &h , &w , S p , & m ) ; 

settextclr(& 1 green , Sgray ) ; 
cx = 18.0 ; cy = 72.0 ; 
movtcurabs ( &cx , Scy ) ; 

text ( "WARNING ! CONTACT") ; 
m = 0 ; 

settext ( &h , &w , &p , &m ) ; 

settextc 1 r (& 1 green , Sblack ) ; 

cx = 35,0 ; cy = 29,0 ; 
movtcurabs ( Sex , Xcy ) ; 

sprintf (s , "%f ", xt ) ; 

text ( s ) ; 

xt = xt + 10.0 / 3600.0 * U * cos(RB) ; 

movtcurabs ( Sex , Scy ) ; 

sprintf ( s , "%f ", xt ) ; 

text(s) ; 

cy = 26.0 ; 

movtcurabs ( Sex , Scy ) ; 

sprintf ( s , "%f " , y t ) ; 

text ( s ) ; 

yt = yt + 10.0 / 3600.0 * U * sin(RB) ; 
movtcurabs( Sex , Scy ) ; 
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spr intf ( s , "%f " , y t ) ; 

text ( s ) ; 

cy = 2 3 . 0 ; 

movtcurabs ( &cx , &cy ) ; 

sprintf (s , "%f " , zt ) ; 

text(s) ; 

mov tcu rabs ( &cx , Scy ) ; 

sprintf (s , "%f " , zt ) ; 

text ( s ) ; 

cy = 20.0 ; 

mov tcu rabs ( Sex , Scy ) ; 

sprintf ( s , "%f ", R ) ; 

text(s) ; 

R = pow(xt*xt + yt*yt + zt*zt , 0.5 ) ; 
movtcurabs ( &cx , &cy ) ; 

spr in tf ( s , " %f " , R ) ; 

text(s) ; 

XD = 50,0 + 45.0 * xt / max_range ; 

YD = 50.0 + 45.0 * yt / max_range ; 

} 

i = getch f ) ; 

cx = 97.0 ; cy = 2.0 ; 
movtcurabs( Xcx , &cy ) ; 

i = getch ( ) ; 

return ; 

* 
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♦♦include "stdio.h'' 
♦♦include "math.h" 
double U , RCFA ; 



MAIN ( ) 

i 

CPA( ) ; 

NUMBER ( ) ; 

exi t ( 0 ) ; 



CPA ( ) 

/* The CPA( ) function computes the CPA between the */ 
/* target and the ship. Two CPA's are computed one */ 
/* is the CPA in X-Y plan ( or horizontal plane ), */ 
/* the other one TCP A ( true CPA ) is talcing the */ 
/* altitude of the target into account. The CPA is */ 
/* converted into TCPA by changing from X— Y plan */ 
/* to Z plan. The parameters needed for this func- */ 
/* tion will be explained when they are declared . */ 



float 


xs , 


ys 


/* 


initial 


position of the 


ship 


*/ 


float 


xt , 


yt 


/* 


initial 


position of the 


target 


*/ 


float 


z t 




/* 


constant 


altitude of the 


target 


*/ 


float 


vs , 


as 


/* 


ve 1 oci ty 


and course of the ship 


*/ 


float 




at 


/* 


ve 1 oci ty 


and course of target 


*/ 


float 


R 




/* 


distance 


between ship and 


1 target*/ 


float 


A 




/* 


angle between horizontal 


and R 


*/ 


float 


xl , 


yl 


/* 


define ship's velocity ve 


ctor 


*/ 


float 


x2 , 


y2 


/* 


define target's velocity 


vector 


*/ 


float 


x , 


y 


/* 


target ' s 


relative velocity vect 


*/ 


float 


B 




/* 


angle of 


x— y vector 




*/ 


float 


C 




/* 


angle between R and x-y v 


ector 


*/ 


float 


CPA 




/* 


CPA in X- 


-Y plan 






/*f 1 o< 


at RCPA 


* 


CPA with 


zt component 




/ 



printf ( "\n PLEASE INPUT X, Y AND Z POSITION OF THE TARGET 
" ) ; 

scanf ( ”%f %f %f " , Sxt , &yt , 8zt ) ; 

printf ("\n PLEASE ENTER ITS SPEED AND COURSE : " ) ; 

scanf ("% 1 f %f " , &U , Sat) ; 

printf ( "\n PLEASE INPUT X AND .Y POSITION OF THE SHIP : " ) 

scanf ( "%f%f " , &xs , &ys) ; 

printf ( "\n PLEASE ENTER ITS SPEED AND COURSE : ; 

scanf ( "%f %f " , Jvs, 8as ) ; 

at = at * PI / 180.0 ; 

printf ("\n AT = %f ", at ) ; 

as = as * PI / 1 80 . 0 ; 
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E = pow ( ( xs-xt ) * ( xs-xt ) + ( ys-y t ) * ( ys-y t ) ,0.5) ; 

printf ( "\n. R = %f " , R ) ; 

A = atan2 ( ( y t-ys ) , ( xt-xs ) ) . ; 

printf ( "\n A = %f " , A) ; 
xl = vs * cos ( as ). ; 

yl = vs * sin(as) ; 

x2 = U * cos ( at ) ; 

y2=U * s i n ( a t ) 

x = xl - x2 + 0.000001 ; /* to avoid division by zero */ 
y = yl - y2 ; 



printf ( "\n X 


= %f 


4-1 

II 


, x , y ) 


B = a tan2 ( y , 


,x ) 


j 




printf ( "\n B 


= %£ " 


, B ) ; 




C = A - B 
printf ( ”\n C 


= %f " 


, C ) ; 





if ( C > 1.57079 I! C < -1.57079 ) 

i 

printf ( "\n NO CPA ! " ) ; 

return ; 

} 

CFA = K * sin(C) ; 

RCPA= pow ( ( CPA * CPA ) + ( zt * zt ) , 0.5 ) ; 

printf ( "\n THE CPA IS : %f ” , CPA ); 
printf ("\n THE TRUE CPA IS : %lf", RCPA ) ; 

return ; 
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NUMBER ( ) 

{ . 

double RI[10] , RT[10] , TA[10] , LA[10] 

double MPHI , PHI , NR , NL , DIST , T , V , RX ; 

double TEMP , RN ; 

double PK [10] , TPK =1.0; 

int N , i ; 



printf ( "\n PLEASE ENTER THE MAXIMUM THEN THE MINIMUM MISSILE 
RANGE"); 

scanf ( " % 1 f % 1 f " , SRX , 2 RN ); 

if ( P.CPA > EX) 
i 

printf ( "\n THE TARGET IS BEYOND THE MISSILE'S RANGE !") ; 

return ; 

} 

printf ("\n ENTER THE SPEED OF MISSILE " ) ; 

scanf ( "%lf " , ScV ) ; 

printf ( "\n ENTER THE DELAY TIME T " ) ; 

scanf ( “%lf " , &T ) ; 



RI [0 J 


= RX ; 


T A [ 0 ] 


= 0 ; 


N 


= 0 ; 


TEMP 


= RCPA/RI [ 0 ] ; 


NPHI 


= asin ( TEMP ) ; 


printf ("\n RCPA = %lf , T = %lf 


DIST 


= T * U ; 


TEMP 


= ( DIST*DIST + RI [0 ]*RI [ 0 


NR 


= pow ( TEMP ,0.5) ; 


TEMP 


= DIST*RCPA/(NR*RI [0] ) ; 


NL 


= asin (TEMP) ; 



printf ( "\n NPHI = %lf 



%lf 
TA [1 ] 
TEMP 
LA [1 ] 
PHI 
RI [1 ] 
RT [1 ] 
N 



DIST 



NR 



. NPHI 
NPHI + NL ; 
U/V*sin(TA[l ] ) 
asin(TEMP) ; 

PI - TA[1 ] - L A [ 1 ] 
RCPA / sin (PHI ) 
U/V*RI [1 ] 

1 



DIST 
, NL 



U = %lf " . RCPA, T, U ) ; 

- 2*DIST*RI [0 ]*cos(NPHI ) ) ; 

= %lf , NR = %lf , NL = 



while ( RI [N ] <= RX ) 
if ( RI [N ] < RN) 

5 

< 

TEMP = RCPA / RN ; 

NPHI = PI - asin ( TEMP ) ; 

TEMP = DIST*DIST + RN*RN - 2*DIST*RN*cos ( NPHI ) ; 

NR = pow ( TEMP ,0.5) ; 

TEMP = DIST*RCPA/(NR*RN) ; 
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ML = asin ( TEMP ) 

TA[N] = NPHI + ML ; 

TEMP = U/V*sin(TA[N] ) ; 

LA [N ] = asin (TEMP) ; 

PHI = PI - TA [N ] - LA[N] 
RI[N] = RCPA. / sin (PHI) 
RT[N] = U/V*RI[N] ; 
if ( RI [N ] >= RX ) 



printf ( "\n NUMBER OF MISSILES : %d - ” , N ) ; 

return ; 

i 

i 

> 

NPHI = PI- - PHI 

= (DIST*DIST + RI[N]*RI[N] - 2*DIST*RI [N ]*cos ( NPHI ) ) 
= pow (TEMP , 0.5 ) ; 



TEMP 

NR 

TEMP 

NL 

N 



= DIST*RCPA/ ( NR*RI [N ] ) 

= asin (TEMP) ; 

= N + 1 ; 

TA [N ] = NPHI + NL ; 

TEMP = U/V*sin(TA[N J ) ; 

LA[N] = asin (TEMP) ; 

PHI = PI - TA [N ] - LA [N ] 
RI[N] = RCPA / sin (PHI) 

RTCN] = U/V*RI [N] ; 

} 

for ( i =0 ; i <N ; i ++ ) 

I 

PK[i ] = 1.233333 - RI[i]/24.0; 
TPK = TPK * ( 1 . 0 - PKCi ] ) 



TPK =1.0- TPK ; 

printf ( "\n RCPA = %lf , RX = %lf , P.N = % 1 f ", RCPA.RX.RN ) 
for ( i =0 ; i <N ; i ++ ) 

printf ( "\nRI [%d]=%lf ,TA[%d] = %lf , PK [%d 3 = %1 f " 

, i , RI [i ] , i , TA [ i ] , i ,PK[i ]) ; 

printf ("\n THE TOTAL PROBABILITY OF KILL IS %lf ", TPK ); 
printf (”\n THE NUMBER OF MISSILES THAT CAN BE FIRED : %d 

N ); 

return ; 
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